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Abstract: We address the nonperturbative structure of topological strings and c = 1 matrix 
models, focusing on understanding the nature of instanton effects alongside with exploring their 
relation to the large-order behavior of the 1/N expansion. We consider the Gaussian, Penner and 
Chern-Simons matrix models, together with their holographic duals, the c = 1 minimal string at 
self-dual radius and topological string theory on the resolved conifold. We employ Borel analysis 
to obtain the exact all-loop multi-instanton corrections to the free energies of the aforementioned 
models, and show that the leading poles in the Borel plane control the large-order behavior 
of perturbation theory. We understand the nonperturbative effects in terms of the Schwinger 
effect and provide a semiclassical picture in terms of eigenvalue tunneling between critical points 
of the multi-sheeted matrix model effective potentials. In particular, we relate instantons to 
Stokes phenomena via a hyperasymptotic analysis, providing a smoothing of the nonperturbative 
ambiguity. Our predictions for the multi-instanton expansions are confirmed within the trans- 
series set-up, which in the double-scaling limit describes nonperturbative corrections to the Toda 
equation. Finally, we provide a spacetime realization of our nonperturbative corrections in terms 
of toric D-brane instantons which, in the double-scaling limit, precisely match D-instanton 
contributions to c = 1 minimal strings. 
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1. Introduction and Summary 



The nonperturbative realm of quantum field and string theories has often been a source of many 
new results and surprises. Another recurrent topic of great interest over the past decades has been 
the large ./V approximation, lately relating gauge and string theories in a nonperturbative fashion. 
Of particular interest to us in this work is the case of the perturbative 1 /N expansion of hermitian 
matrix models, whose nonperturbative corrections are exponentially suppressed as exp {—N). In 
the double-scaling limit these models describe noncritical or minimal (super)string theories, and 
the nonperturbative structure of the matrix model is related to that of the corresponding string 
theory Q: the exp (— N) contributions are instanton effects in the matrix model j|, [| and 
they are interpreted as D-brane configurations in the string theoretic description Eh q|. Of 
course the study of matrix models is not confined to the vicinity of their critical points and 
one may also study nonperturbative effects away from the double-scaling limit. The interesting 
point is that off-critical matrix models may be dual to topological string theories. For instance, 
this happens in the case first suggested by 0, where some off-critical matrix models describe 
the topological string B-model on certain non-compact Calabi-Yau (CY) backgrounds, with 
the string genus expansion (in powers of the string coupling, g s ) being identified with the 1/N 
matrix model expansion; and it is also the case for topological strings with a Chern-Simons 
dual, first studied in [||, (9|. This turns out to be a more general statement, as it was later 
shown in ]l0| , 11, |l^] that topological string theory on mirrors of toric manifolds also enjoys a 



dual holographic description in terms of off-critical matrix models. It is thus evident that fully 
understanding the nonperturbative structure of matrix models, both at and off criticality, will 
have many applications in both minimal and topological string theories. 

Recently [13, 14, O] there has been significative progress in understanding and in quanti- 



tatively computing nonperturbative effects in matrix models away from criticality. In 13], and 
building upon double-scaled results g |3| |l6|, [l7|], off-critical saddle-point techniques were de- 
veloped in order to compute instanton amplitudes (up to two loops) in terms of spectral curve 
geometrical data. This work focused upon one-instanton contributions in one-cut models, and 
in [15 1 an extension to multi-instanton contributions, again in one-cut models, was obtained, 



starting from a two-cut analysis. Extensive checks of the nonperturbative proposals in these 
papers were also performed, by matching against the large-order behavior of the 1/N expansion. 



Another approach to multi-instanton amplitudes was developed in [14], this time around based 



on orthogonal polynomial methods, via the use of trans-series solutions to the string equations 



[18]. Further progress along these lines recently led to the proposal of [19], where a proper 
nonperturbative definition of a modular-invariant holomorphic partition function was presented, 
which was also shown to be manifestly background independent. Remarkably, many of the re- 



sults uncovered in [O, 14, [Lq] appear to extend beyond the context of matrix models; e.g., in 



cases where the theory is controlled by a finite-difference equation — such as the string equation 



[18] for matrix models — it is possible to compute nonperturbative effects and relate them to the 



large-order behavior of the theory. This is the case of Hurwitz theory [13], which is controlled by 
a Toda-like equation, and also the case of topological strings on the background considered in 
pOjl 1 . However, all models considered in the aforementioned articles lie in the universality class 
of two-dimensional gravity, with c = 0, and methods that have been worked out in this case can- 
not be applied in a straightforward fashion to the case of topological strings in the universality 
class of c = 1. In view of this, it is necessary to develop new techniques in order to approach 



1 A local CY threefold given by a bundle over a two-sphere, X p = 0(p — 2) © C(— p) 4 P 1 , p 6 Z, which may 
be regarded as a quantum group deformation of Hurwitz theory; see feci] for further details. 
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nonperturbative effects in models which belong to the universality class of the c = 1 string at 
the self-dual radius. 

Let us be a bit more specific about the nature of the string perturbative expansion and the 
type of nonperturbative contributions we shall be looking for. Topological strings, much like 
physical string theory, are perturbatively defined in terms of two couplings, ol and g s , as 2 

+00 

F(g s ;{t t }) = ^9 2 s 9 - 2 F g (t l ), (1.1) 

where F = log Z is the free energy and Z the partition function, and where the fixed genus free 
energies F g (ti) are themselves perturbatively expanded in a'. In some sense the a' expansion 
is the milder one: it has finite convergence radius, with this radius given by the critical value 
of the Kahler parameters where one reaches a conifold point in moduli space. As it turns out, 
the problem of finding a nonperturbative formulation of the A-model free energy, in a' , may be 
reduced to that of solving the mirror B-model description, where topological string amplitudes 
become exact in a' . In this way, the A-model solution is found by translating B-model amplitudes 
back to the A-model, by means of the mirror map. This topic has been extensively studied in 



the literature and we refer the reader to the recent developments [11, 12] and references therein. 

The situation gets more complicated as one tries to go beyond perturbation theory in g s . 
In this case, one is immediately faced with the familiar string theoretic large-order behavior 



F g ~ (2<?)! rendering (1.1) as an asymptotic expansion H|. In this case, one expects nonper- 
turbative corrections of order ~ exp(— l/g s ), and an adequate nonperturbative formulation of 
the theory must encode all these corrections. As described above, there are certain cases — such 
as the backgrounds considered by Dijkgraaf and Vafa [j?], or models with a dual Chern-Simons 
interpretation — where topological strings have a holographic matrix model description, with the 
matrix model large N expansion reproducing the topological string genus expansion. In these set 
of backgrounds one would be tempted to use the finite N matrix model free energy as the g s non- 
perturbative definition of topological string theory 3 . In order to establish this result, one must 
first understand how the finite N matrix model would encompass all nonperturbative contribu- 
tions ~ exp (— l/gs)- This situation is clear for minimal strings, realized in the double-scaling 
limit of hermitian matrix models: the nonperturbative effects associated to the asymptotic nature 
of the genus expansion are implemented via eigenvalue tunneling effects in the dual matrix model, 
and are interpreted in the continuum formulation in terms of Liouville branes in spacetime jl], || . 
For topological strings, a similar understanding has been achieved in the case of the local curve 
flO| , [DJ, where a matrix model description is available |{2C[| . In this case, the nonperturbative 
effects associated to the asymptotic behavior, or large-order behavior, have again been matched 
to instantons arising from matrix eigenvalue tunneling, and a spacetime interpretation in terms 
of domain walls has been provided |l3[| . 

However, there are several cases where this paradigm seems not to apply, at least not in 
a straightforward fashion. It is our goal to address such issues in the present work in the 
prototypical example of the resolved conifold, but also encompassing matrix models in the c = 1 
universality class. Topological strings on the resolved conifold are holographically described 
by the Chern-Simons matrix model, but there are now no obvious instantons associated to 
eigenvalue tunneling as the Chern-Simons potential has no local maxima outside of the cut, 
where the eigenvalue instantons could tunnel to. This problem, which was not an issue in any of 

2 Recall that in the A-model the {ti} are Kahler parameters while in the B-model they are complex parameters. 
3 Other nonperturbative completions, provided by a holographic dual, have been proposed in pi] ]. 
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the previously mentioned examples, also appears in other matrix models, such as the Gaussian 
and Penner models; all of them in the c = 1 universality class. One may then ask where do 
nonperturbative corrections arise from, or what exactly controls the large-order behavior of the 
1/N perturbative expansion in these models. We shall answer these questions in this paper. 
One way out is to directly compute the (would-be) instanton action that controls the large- 
order behavior of the perturbative expansion, by means of a standard Borel analysis (see, e.g., 
||22|| ). At first this may look like a formidable task, as one may expect the topological string 
genus expansion to be rather complicated, not amenable to a Borel transform. Happily, the free 
energies of all cases we consider enjoy a Gopakumar-Vafa (GV) integral representation |23|, |24j 
which allows for an exact location of the singularities in the Borel complex plane controlling 
the divergence of the asymptotic perturbative series, i.e., the instanton action This is the 
topological string generalization of a celebrated c = 1 string result [25|. Interestingly enough, 
this integral representation may also be regarded as an one-loop Schwinger integral f24j| , thus 
providing a spacetime interpretation of these nonperturbative effects; as already pointed out in 
[24] they control the pair-production rate of BPS bound states. As we shall later see, these 
results — which we further identify as Stokes phenomena of the finite N partition function — will 
also allow us to explain the nonperturbative contributions as one-eigenvalue effects in the matrix 
model picture. We find, from a saddle-point analysis, that the c = 1 nonperturbative effects arise 
due to the multi-valued structure of the effective potential (as preliminarily suggested in (^6|); a 
different picture from that of matrix models in the universality class of 2d gravity plus matter, 
where the one-eigenvalue tunneling occurs from a metastable minimum to the most stable one. 

This paper is organized as follows. We begin in section 2 by reviewing the main ideas behind 
our subsequent work. This includes the definition of the Borel transform and the relation between 
instantons and the large-order behavior of perturbation theory, both related to the existence of 
a nonperturbative ambiguity in the calculation of the free energy. In this section we also discuss 
the Schwinger effect, where one actually has a physical prescription to define the inverse Borel 
transform, which will turn out to be the case for topological strings and c = 1 matrix models via 
the GV integral representation of the topological string free energy. In section 3 we then move on 
to presenting the matrix models we shall be focusing upon. We review some of their properties, 
such as their spectral curves and their perturbative genus expansions, and also obtain expressions 
for their exact, finite N partition functions and holomorphic effective potentials, both of which 
play important roles in sections to come. In this section we also discuss the double-scaling 
limit of these models and show how they relate to FZZT branes. Section 4 presents one of the 
main topics in this paper, the Borel analysis of the Gaussian, Penner and Chern-Simons matrix 
models. We show how to obtain Schwinger-like integral representations of the free energy, via 
Borel resummation, and how the correct identification of the leading poles in the complex Borel 
plane leads to the one-instanton action in all our examples. We further show in this section that 
while the all-loop multi-instanton amplitudes precisely reconstruct the perturbative series, the 
one-instanton results control the large-order behavior of perturbation theory. We then move on 
to another of our main topics in section 5, namely the issue of Stokes phenomena. We recall 
how to obtain Stokes phenomena for integrals with saddles via hyperasymptotic analysis, and 
perform a detailed calculation for the Gamma function. This extends to the Barnes function 
and, in this way, allows us to identify instantons with Stokes phenomena as we reproduce the 
results we have previously found in section 4, out of hyperasymptotic analysis. In section 6 
we provide a semiclassical interpretation of our instantons via eigenvalue tunneling, where this 
tunneling is now associated to the existence of a branched multi-sheeted structure in the relevant 
holomorphic effective potentials. Indeed, simple monodromy calculations reproduce our results 
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for the multi-instanton action straight out of this interpretation. We further show in this section 
how to interpret our instantons in spacetime, from the point of view of ZZ branes. In section 
7, we discuss the trans-series approach to c = 1 matrix models and how it further validates our 
results. Finally, we conclude in section 8 with an outlook and future prospects. We also include 
two appendices, one dedicated to the study of the monodromy structure of the polylogarithm, 
and the other dedicated to the Cauchy dispersion relation, in the case of more general topological 
string theories than the ones we address in this paper. 



2. Asymptotic Series, Large Order and Topological Strings 

We start by reviewing some useful facts concerning asymptotic series, the relation of their large- 
order behavior to nonperturbative effects, as described by instantons or by the Schwinger effect, 
and put them in the context of topological string theory as we wish to study in the present work. 
For an introduction to these topics with applications in quantum mechanics and quantum field 
theory, we refer the reader to |22j and references therein. 

Let us consider the perturbative expansion of some function, F(z), with z the specific per- 
turbative expansion parameter, 

F(z)~^F n z™. (2.1) 

n=0 

In many interesting examples one may infer that, at large n, the coefficients behave as F n ~ (/?n)!, 
thus rendering the series divergent. As an approximation to the function F(z), the asymptotic 
series ( p.ip must necessarily be truncated. As such, one is faced with an obvious problem: how to 
deal with the fact that the perturbative expansion has zero convergence radius? In particular, if 
we do not know the function F(z), but only its asymptotic series expansion, how do we associate 
a value to the divergent sum? The best framework to address issues related to asymptotic series 



is Borel analysis. One starts by introducing the Borel transform of the asymptotic series (2.1) as 

+oo 

m(o = Y,^ n > (2-2) 

which removes the divergent part of the coefficients F n and renders 0[.F](£) with finite conver- 
gence radius. In particular, if F{z) originally had a finite radius of convergence {i.e., if it was 
not an asymptotic series), B[F](^) would be an entire function in the Borel complex £-plane. In 
general, however, B[F](^) will have singularities and it is crucial to locate them in the complex 
plane. The reason for this is simple to understand: if B[F](^) has no singularities for real positive 
£ one may analytically continue this function on M + and thus define the inverse Borel transform 
by means of a Laplace transform as 4 

F(z) = / dsB[F](z S P)e- s . (2.3) 
Jo 

The function F(z) has, by construction, the same asymptotic expansion as F(z) and may thus 
provide a solution to our original question; it associates a value to the divergent sum (^]l|). If, 
however, the function S[i ? ](£) has poles or branch cuts on the real axis, things get a bit more 
subtle: in order to perform the integral (|2.3D one needs to choose a contour which avoids such 



'For simplicity, we are assuming z £ R + in this expression. 
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singularities. This choice of contour naturally introduces an ambiguity (as we shall see next, 
a nonperturbative ambiguity) in the reconstruction of the original function, which renders F(z) 
non-Borel summable 5 . As it turns out, different integration paths produce functions with the 
same asymptotic behavior, but differing by exponentially suppressed terms. For instance, in the 
presence of a singularity at a distance A from the origin, on the real axis, one may define the 
integral ( |2.3| ) on contours C±, either avoiding the singularity from above, and leading to F+(z), or 
from below, and leading to F^(z). One finds that these two functions differ by a nonperturbative 



term 22 



F+(z)-F-(z) ~ ie (2.4) 

In certain cases, e.g., when one has a Schwinger representation for the function [27, 28], there is a 
natural and rigorous way to define the integral (^1|) on a contour which avoids the singularities, 
and which also allows for a physical interpretation of the nonperturbative contributions. 

So far our discussion has been rather general. However, it takes no effort to figure out the 
physical relevance of our discussion: divergent series are almost ubiquitous in physics and appear 
basically each time we approach an interesting problem in perturbation theory p2]|. A typical and 



extensively studied case in quantum mechanics is the anharmonic oscillator (see, e.g., |£9, 3C , P2|). 
Herein, the ground state energy may be computed in perturbation theory — as a power series in 
the quartic coupling — and one finds that it is analytic in all the (coupling constant) complex plane 
except for a branch cut on the negative real axis, associated to the instability of the potential 
which becomes unbounded for negative values of the coupling. This instability is reflected by the 
fact that the series is, as expected, asymptotic. In particular, one can perform a Borel analysis as 
above and discover that the Borel transform of the ground state energy has singularities on the 
positive real axis, leading to an ambiguity of order ~ ie~ 1//fil , with g the quartic coupling constant. 
In this simple quantum mechanical example the nonperturbative ambiguity has a clear physical 
interpretation: it signals the presence — at negative g — of instantons mediating the decay from 
the unstable to the true vacuum, via tunneling under the local maximum of the potential. 

What these ideas illustrate is that by means of a purely perturbative analysis, i.e., finding the 
singularities of the Borel transform of the original perturbative series, it is possible to learn about 
nonperturbative effects — at least the intensity of the nonperturbative ambiguity (but we shall say 
more on this in the following). In some examples it is possible to independently compute these 
nonperturbative terms directly, e.g., using WKB methods or computing the path integral around 



non-trivial (subdominant) saddle points [22|. In these examples one may then proceed in the 
opposite way from above and obtain information on the large-order behavior of the perturbative 
expansion out of the nonperturbative data. This is what we shall illustrate next. 

2.1 From Instantons to Large Order and Back 

In physical applications, the factor A appearing in ( |2.4D is the one-instanton action (see, e.g., 
[22]). Let us make this relation between instantons and the large-order behavior of perturbation 



theory a bit more precise, as it will play a crucial role in our later analysis. Consider a quantum 
system whose free energy is expressed as a perturbative expansion in g, the coupling constant 6 , 

^ 0) (<7) = £/f<A (2.5) 

k=0 

5 Strictly speaking, the function is said not to be Borel summable if different integration contours yield different 
results. It may still be the case that, in spite of having singularities in the real axis, all alternative integration 
contours yield the same result. 

6 Here and below the index (£) labels the £-instanton sector, so that (0) labels the perturbative expansion. 



-6- 



The series fl2.5|) will generically be asymptotic, with zero radius of convergence. This is naturally 
associated to a branch cut of F(g) in the complex g-plane, located in the negative real axis 
and associated to instanton effects (just like in the anharmonic oscillator example above). The 
function F(g) is expected to be analytic otherwise. In fact this is saying that our quantum system 
should actually be thought of as an asymptotic formal power series in two expansion parameters, 
g and exp (--], see [14] for a discussion in the matrix model context. The appropriate expansion 



of the free energy is thus [14] 



, oo 

IA 



F(g) = ^C*F<%), F(%) = (2.6) 



Q 



Here, C is a parameter corresponding to the nonperturbative ambiguity. Also, A is the one- 

(£) 

instanton action, b a characteristic exponent and ft is the fc-loop contribution around the 



-instanton configuration. Typically, the coefficients are factorially divergent for any £ [22], 



in which case we may think about the (.£+1)— instanton sector as the nonperturbative contribution 
related to the asymptotic nature of the loop expansion around the ^-instanton sector. 

A standard procedure then relates the coefficients of the perturbative expansion around the 
zero-instanton sector, fj^\ with the one-instanton free energy as follows. The discontinuity of 
the free energy across the branch cut (associated to the instability of the theory for negative g) 
is expressed, at first order, in terms of the leading instanton expansion ( |2.6| ) 

Disc F(g) = lim F(g + ie) - F(g - ie) = 2i Im F(g) = F (1) (g) + ■■■ . (2.7) 
At the same time, we may use the the Cauchy formula to write 



2vri 7_ oc w-g / (oo) 2vri w 

In certain situations, e.g., in the aforementioned anharmonic oscillator example |29( [, it is possible 
to show by scaling arguments that the last integral in the expression above does not contribute. 
In such cases, (2.8) provides a remarkable connection between perturbative and nonperturbative 



expansions. Using the perturbative expansion (2.5) and the leading one-instanton contribution to 



the discontinuity Disc F(g) ~ F^ , one may obtain from the Cauchy formula ( |2.8] ) the following 
large order (or large k) relation 

(0) [ +ao dz_ FW(z) T(k + b) ^ T( k + b-n) (1) 
Jk J 2vri ~ 2irA k + b ^ V (k + b) In+1 ' [ ' 

n=0 

This explicitly shows that the computation of the one-loop one-instanton partition function 
determines the leading order of the asymptotic expansion for the perturbative coefficients of the 
zero-instanton partition function. Higher loop corrections then yield the successive r corrections. 
Furthermore, instanton corrections with action A 1 > A, where we have in mind multi-instanton 
corrections with action £A, I > 2, will yield corrections to the asymptotics of the coefficients 
which are exponentially suppressed in k. 

For the cases we shall consider in this work, namely matrix models and string theory, one 



finds genus expansions as in (1.1), with Fg ^ ~ (2<?)!, so that the relation (^j|) gets slightly 
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re-written as follows (see, e.g., |13|]). Begin with the free energy in the zero-instanton sector, 
g^F^ (g s ). Setting z = g^, the one-instanton path integral then yields a series of the form 

+00 

zF^(z) = -is-e-v* $>**&\(t). (2.10) 
Following a procedure analogous to the one above, where one further assumes that the stan- 



dard dispersion relation (2.8) still holds, it follows for the zero-instanton sector perturbative 
coefficients 



F (0 )m _ f + °° dz zFW (z) r (2g + b) +^ r (2g + 6 - fe) (1) ft 

Again, the computation of the one-loop one-instanton free energy determines the leading order 
of the asymptotic expansion for the perturbative coefficients of the zero-instanton free energy. 
Higher loop corrections then yield the successive corrections. One should further notice that 
recently, in [jl3| , [lj, [ljj, the relation ( [2.11D has been tested in several models and rather conclusive 
numerical checks have confirmed its validity. 



2.2 The Schwinger Effect and a Semiclassical Interpretation 

As discussed above, a nonperturbative ambiguity — typically associated to instantons, from a 
physical point of view — can arise when defining the integration contour for the inverse Borel 
transform. However, this is not always the case: we shall now review an example where a 
prescription to define the inverse Borel transform naturally arises, together with a physical in- 
terpretation for the nonperturbative contributions p8j. This is the Schwinger effect [27]. 



The one-loop effective Lagrangian describing a charged scalar particle, of charge e and mass 
m, in a constant electric field, E > 0, has an integral representation given by ]27j (see, e.g., [31] 
for a recent review) 



e 2 E 2 f +oo ds / x 1 s 



■2 



which admits the weak coupling expansion 

4 +00 , 9 p \ 2n+4 



C^JH-y { - ir _^±i_ [™\ . (2.13) 

167r n^o (2ra + 4)(2n + 3)(2n + 2) \ m 2 / 



In here we used the shorthand B<m = 1 2 ^-i -^2m with i?2n the Bernoulli numbers. Since one 
may further relate the Bernoulli numbers to the Riemann zeta function via 

^ = (-D (B+1, ^(2n)I ( 2 - 14 ) 

it becomes evident that Bm diverges factorially fast. In this case, if one first writes the expansion 

as 

+00 

£~5> 2n+4 x 2 " +4 , (2.15) 

n=0 
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with x = ^jr, it follows that at large n one has 02^+4 ~ (2ra + 1)! rendering this perturbative ex- 
pansion asymptotic — and actually non-Borel summable as we shall see next. Indeed, computing 
the Borel transform it follows 

from where one immediately notices that the Schwinger integral representation of the effective 
Lagrangian ( [2.12 ) is essentially the inverse Borel transform 



£(,) = y o ^B[£](xt)e* = I ^y o _^_---_je*-. (2.17) 

Of course so far we still have a nonperturbative ambiguity to deal with: in order to perform the 
integration on the real axis one still needs to specify a prescription in order to avoid the poles 
at s = nir, n G N. This introduces the usual ambiguities leading to exponentially suppressed 
contributions to the effective Lagrangian. The novelty in this case is that there is now a natural 
way to address the integration avoiding the singularities in an unambiguous way psf . 

As it turns out, the contour of integration needs to be deformed in such a way that the 
integral picks up the contributions of all the poles as if the real axis is approached from above, 
tantamount to a +ie prescription; and this is the requirement which is dictated by unitarity [28]. 



As such, one has a physical principle behind the unambiguous choice of contour. Furthermore, 
the Lagrangian develops an imaginary part which is simple to compute by summing residues 7 , 
and which cannot be seen to any finite order in perturbation theory, 



, 1 (eE\ 2 +£?(-l)— 1 / m'\ , s 

■^=^(v) e4-"Hf ■ (2 - 18) 

n=l 



an expression with an evident multi-instanton flavor [32|, as in ( |2.6| ). Besides the appropriate, 
physical prescription to perform the integration, and as such unambiguously compute the non- 
perturbative contributions to the Lagrangian, the Schwinger effect gives us something else: a 
physical interpretation of this imaginary part. Indeed, the imaginary part of the effective La- 
grangian ( p. 18 ) is precisely the pair-production rate, or probability per unit volume for pair 



creation, for scalar electrodynamics in a constant electric field [27]. In other words, the above 
unitary +ie prescription for the integration contour guarantees that this probability is a positive 
number between zero and one (which basically demands ( 2.18] ) to be real and positive). 



Another interesting illustration of the Schwinger effect, which will be of particular relevance 
in our subsequent discussion on topological strings and matrix models, is the case of a constant 



(euclidean) self-dual electromagnetic background |33], |3J, [35[ , satisfying 



F/jiu — ~*F[xv — 2 ^/w / pA k ^ • (2.19) 

Following 1 31], we introduce T 2 = \F llv F iiV and the natural dimensionless parameter 7 = . 
In this case, the one-loop effective Lagrangian describing a charged scalar particle is now given 

by 11,0 

e 2j-2 H-oo ds / 1 1 1\ _2s 



7 Notice that the residue at s — precisely vanishes. 
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admitting the weak coupling expansion 



4 +°° d / n t\ 2n+2 

C — V 2n+2 ( I (121) 

16vr 2 ^ 2n (2n + 2) V m 2 / ' 1 " 1 

n=l ' 

Notice that there are two possible self-dual backgrounds [33, 35]: a magnetic-like background 



with T real, in which case ( 2.2 1]) h as an alternating sign; and an electric-like background with 



T imaginary, in which case (|2.21 ) is not alternating. If one carries through a Borel analysis 
similar to the previous one, where we studied the case of constant electric field, one may further 
notice that while the alternating (magnetic) series has Borel poles on the positive imaginary 
axis, the non-alternating (electric) series has the Borel poles on the positive real axis, making 
both situations rather distinct on what respects evaluating the inverse Borel transforms. For 
the non-alternating (electric) series it is also possible to see that the aforementioned unitarity 
prescription will pick a Borel contour leading to a nonperturbative imaginary contribution to the 
Lagrangian, unambiguously given by 

em 2 J 7 r-^ / 2-7T 7 \ _2™ 

Im£ = _yw_ + _L]e y. (2.22) 

32vr 3 n 2 ' 



This expression can similarly be obtained by first considering the magnetic series, reflecting the 
integrand in (2.20) to the negative real axis in order to obtain an integral over the entire real 



line, and then deforming this integration contour such that it just incloses all the poles in the 
positive imaginary axis (a contour surrounding ilR + ). The resulting integral will then produce a 
sum over residues which, upon "Wick rotation" of the dimensionless coupling 7 — > ry, leads to 
the same expression as above, ( p. 22 ). As we shall see in the course of this work, this expression is 



also at the basis of the nonperturbative structure of topological strings and c = 1 matrix models. 

Another important feature of the Schwinger effect, that we shall further explore later on, is 
the fact that in the presence of a constant electric field the pair-production process can be given a 
semiclassical interpretation in terms of a tunneling process, where electrons of negative energy are 
extracted from the Dirac background by the application of the external field [pq| . The motion 
under the potential barrier, classically forbidden, is considered for imaginary values of time, 
allowing for a computation of the tunneling probability corresponding to the pair-production 
rate as 

u^e" 21 ™ 5 , (2.23) 

where ImS is the imaginary part of the action developed during motion under the barrier. In 
here, a crucial point is that a particle in a sub-barrier trajectory satisfies the classical equations 
of motion. One may then use standard classical mechanics of a relativistic particle in order to 
describe this process. Indeed, energy conservation 



£ = ±^p 2 + m 2 - eEx, (2.24) 

together with the equation of motion dtp = eE, allow for an immediate re-writing of the action 
(for 5 = 0) as: 

2 

S{p) = tt~^ \/p 2 + m 2 - —— log (p + \J p 2 + m 2 ) . (2.25) 
let, leh V / 

Notice that, because of the logarithm, the action is a multi-valued function. 

The spectrum of possible values for the energy is displayed in Figure [l|. A potential barrier 

separates the lower continuum of negative-energy states (the minus sign of y/p 2 + m 2 ) from the 
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Rep 



Figure 1: On the left, the energy as a function of the position. The blue (orange) region corresponds to 
the lower (upper) continuum of particles with negative (positive) energy. The white region is classically 
forbidden. On the right, the variation of the imaginary momentum in the sub-barrier motion. 



upper continuum of positive-energy states (the plus sign of \/p 2 + m 2 ). Sub-barrier motion 
between points A\ and A2 will start at A\, where t = = p, and corresponds to the variation 
of the imaginary time/momentum along the path A1BA2, while the real part of the energy 
remains constant. Indeed, at the classical turning point B we will have t = im/eE and p = im, 
which corresponds to a square-root branch point of the function S(p). The motion ends back at 
t = = p in point A2, as shown in Figure |l|. In this case, we see that the sub-barrier trajectory 
correspond to an increment of the imaginary part of the action as 

A AlB A 2 S = ImS, (2.26) 

which we may compute as the shift of the multi-valued function S(p) as we move in-between the 
sheets of the logarithm. In fact, it is rather simple to realize that the value of S(p) on a generic 
sheet differs from its value on the principal sheet, S*(p), by 

2 

77-777 

S( p ) = S*(p) + in^ir, nGZ. (2.27) 
eh 

In the illustration above we went (once) "half way" around the branch cut [—im, +im] in which 
case the shift in the action is given by 

2 

7TTYI 

A AiBA 2 S=—. (2.28) 

For a generic sub-barrier motion, corresponding to a repeated wandering of the particle between 
the turning points A\ and A2 we may write 

r 2 
2ImS=<t yV +m 2 = n TIL n G Z, (2.29) 

hn eE 

where j n is a contour encircling n-times the branch cut of the action. One may notice [p^] that 



this result is in complete agreement with the Schwinger computation result ( 2.18 ). 

Naturally, this semiclassical argument may be refined in order to reproduce the pre-factors 
of the exponential term, and also so as to include the effect of a magnetic field. The point we 
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wanted to make is that, from a semiclassical perspective, the instanton action describing the 
pair-production rate as a tunneling process may be computed via the branch cut discontinuities 
of the multi-valued function S(p). This is a technique that we shall deploy later on in order to 
provide for a semiclassical interpretation of nonperturbative effects in c = 1 matrix models in 
terms of eigenvalue tunneling. 

In this paper we shall apply the techniques we have just described, Borel analysis, instanton 
calculus, and the Schwinger effect, in order to study the nonperturbative structure of topological 
strings and c = 1 matrix models. As such, we now turn to topological string theory with emphasis 
towards the integral representation of its free energy. 

2.3 The Topological String Free Energy 

Asymptotic series and the Schwinger integral representation also appear in the context of topo- 



logical string theory (see, e.g., |]37[ for an introduction). Let us start by describing the free energy 
of the A-model. The closed string sector of the A-model is a theory of maps <ft : T, g — > X from 
a genus-g Riemann surface, S g , into a CY threefold X, which may be topologically classified by 

their homology class /? = [0*(E g )] G H2(X,Z,). One may expand /3 = X^li^ n * on a 
basis [Si] of H2(X,Z), with associated complexified Kahler parameters tj. 

The topological string free energy has a standard genus expansion in powers of the string 
coupling g s , as in (|1.1[) , which in the large-radius phase (i.e., for large values of the Kahler 



parameters, in units of a') becomes 

F(g s ;{ti}) = Y,9 2 s 9 ' 2 F g (ti), F g (U) = J2 N ^Q - ( 2 -30) 

9=0 /3>0 

Here, the sum over (3 is a sum over topological sectors or, equivalently, over world-sheet instan- 
tons. We have further introduced Qi = e~ t{ , with Q 13 denoting HiQ?'' an( ^ we have chosen 
units in which a' = 2tt. The coefficients N g p are the Gromov-Witten invariants of X, counting 
world-sheet instantons, i.e., the number of curves of genus g in the two-homology class j3. 

The expansion in world-sheet instantons in ( |2.30| ), regarded as a power series in , gener- 
ically has a finite convergence radius, t c , that can be estimated from the asymptotic large j3 



behavior of Gromov-Witten invariants [38] 



N g ,p ~ ( Q(7-2)(l-9)-l e ^c ) p _^ +oc . (2.31) 

In here, 7 is a critical exponent. At the critical value of the Kahler parameter, t c , the so-called 
conifold point, the geometric interpretation of the A-model large-radius phase breaks down 
and the topological string free energy undergoes a phase transition to a non-geometric phase, 
nonperturbative in a'. One may characterize the theory by its critical behavior at the conifold 
point. In particular, one can consider the following double-scaling limit 

t->t c , g s ^0, [i= fixed. (2.32) 

9s 

In this case, the double-scaled free energy is universal, as first noticed in and reads 

1 1 + °° B 

FbsM = Fc=i((-i) = ^ 2 log^ - - log/, + 2 g(2g 9 -2f 2 ~ 2g ' ^ 
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where F c= \((j) is the all-genus free energy of the c = 1 string at the self-dual radius (for a 
review on these issues see, e.g., |4Q[|). The critical behavior ( |2.33 ) has been checked in many 



examples, such as [41, 42]. Furthermore, in |2Ql , it has been shown that certain local CYs have a 
critical behavior which is in the universality class of 2d quantum gravity, i.e., they have 1 = —\- 
Another feature to notice is that the above genus expansion depends on the alternating Bernoulli 
numbers and, thus, is alternating for real \x. 

Of particular interest to our present work is the fact that the free energy F c= i(fi) has a 



Schwinger-like nonperturbative integral formulation 25, 43], given by 



which coincides, after an appropriate identification of the parameters, with the one-loop effec- 
tive Lagrangian for a charged particle in a constant self-dual background, ( 2.2C| ) . This means 



that F c= i{ji) enjoys an asymptotic weak coupling expansion as in ( 2.2l| ) and further develops a 



nonperturbative imaginary contribution akin to (2.22). In this line of thought, the exploration 



of Schwinger-like integral representations for the free energies of topological strings and c = 1 
matrix models is one of the main topics in this paper. 

2.4 A Schwinger Gopakumar Vafa Integral Representation 

As should be clear by now, Schwinger-like integral representations for the free energy are bound 
to play a critical role in our analysis. Happily, for topological string theory, such representations 



have been provided by Gopakumar and Vafa in f23| , 24, [44|. These works explored both the 
connection of topological strings to the physical IIA string, as well as the duality between type 
IIA compactified on a CY threefold, at strong coupling, and M-theory compactified on the 
same CY times a circle, in order to relate topological string amplitudes to the BPS structure of 
wrapped M2-branes and thus re-write the topological string free energy in terms of an integral 



representation. The final result in [44] for the all-genus topological string free energy, on a CY 
threefold X , is 

I-+CO J , a\2r — 2 f 2-7TS \ 

F x (g s ;{U})= £ nW(X) -( 2sin o) exp (d-t + im)). (2.35) 

{di},r,m 

Let us explain the diverse quantities in this expression. The integers (X) are the GV 
invariants of the threefold X. They depend on the Kahler class di and on a spin label r. Later 
on we shall be focusing on the case where X is the resolved conifold, for which there is only one 
non-vanishing integer, = 1. The combination Z = d ■ t + \ m represents the central charge 
of four-dimensional BPS states obtained in the following fashion [44]. Start with M-theory 
compactified on X x S 1 and consider the BPS spectrum of M2-branes wrapped on cycles of the 
CY threefold with fixed central charge A = d ■ t = yj^im-f^OY.Z) diti, with di as above and ti the 
complexified Kahler parameters. The mass of the wrapped M2-branes is 2ttA. Upon reduction on 
S 1 each BPS state may have in addition an arbitrary (quantized) momentum m around the circle, 
leading to BPS states of central charge Z and mass 2ttZ. Notice that these four-dimensional 
BPS states contributing to the topological string free energy may be understood, from a IIA 
point of view, as bound states of D2 and DO-branes, and it is the physics of this system which 
can be related to a Schwinger-type computation and thus to the above integral representation 
(in fact, thanks to the N = 2 supersymmetry in the problem, the Schwinger calculation one 
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has to perform in this context turns out to be equivalent to that of a vacuum amplitude for a 
charged scalar field in the presence of a self-dual electromagnetic field strength, as in ( |2.20| )). 
Furthermore, the integer m, associated to the winding around S , counts the number of DO 
branes in the D2D0 BPS bound state. This should make ( 2.35j ) clear. 



One may also recover the perturbative genus expansion from this integral representation. 
Using the familiar identity 



ex P ( -27rim — J = 5 ( n ) , (2.36) 



with 5(x) the Dirac delta function, one may explicitly evaluate the sum over m in ( 2.35 ) and 
thus obtain, after the trivial integration over s, 

+oo +oo +oo 1 _ „ 

^( 5s ;{ta) = EE^ ) wE^( 2sin! | i ) " ( 2 - 37 ) 

r=0 di=l n=l 

This result expresses the topological string free energy, on a CY threefold X, in terms of the GV 



integer invariants p3|, 24, 04j. To be completely precise, it is important to notice that in order 



to recover the full topological string free energy one still has to add to (2.37) the (alternating) 
constant map contribution ^(| 

* ta , = g^„(*)J=£!!^, (2.38, 

where Xk(X) = 2 (h 1 ' 1 — /i 2 ' 1 ) is the Euler characteristic of X. This term can also be given a 
Schwinger-like integral representation. From the point of view of the duality between type IIA 
and M-theory, this amounts to considering only the contribution arising from the DO-branes, or 



Kaluza-Klein modes. The result is [24| 

"+0O 



Fk( 9s ) = Ixk(X) W £ e— L = \ XK {X) Y, ~ 

mfZ \<4/ n=l V ^ 



(2.39) 



In this paper we shall mainly consider the resolved conifold, a toric CY threefold for which 
&\mH2(X , Z) = 1 and thus the only non-vanishing integer GV invariant is = 1. In this case, 
the GV integral representation ( [2.35D immediately yields 



an expression which carries a Schwinger flavor, as we have seen above. It is also important 
to point out that the case of r = is the only one in which the integrand of the GV integral 
representation will have "interesting" poles, i.e., poles of the sine function on the real axis. When 
r > the only poles of the integrand will be at zero and oo in the Borel complex plane. So, in 
particular, when studying more complicated CY threefolds where there is a sum over r > 0, it 
will always be the contribution from GV invariants with r = which will be the most relevant 
for the Schwinger analysis we shall carry through later in the paper and, as such, the case of the 
resolved conifold is a prototypical example for those situations. From the previous expression it 
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is also simple to obtain the perturbative expansion, by summing over m as previously described, 
and one obtains 

, +00 

n=l \ z J 

By carrying through this sum, expanding in powers of g s , and adding the constant map contri- 
bution, one finally obtains the resolved conifold genus expansion as 

F g (t) = (- 1 ) 3 I^A- 2 | \B 2g \ t) 

9K) 2g{2g-2)(2g-2)\ ^ 2g{2g-2)\ 3 29 { 1 1 } 

with Lip(x) the polylogarithm function. We shall later see how a Borel analysis allows for a 
nonperturbative completion of this expansion and moreover how to relate this nonperturbative 
completion to the large-order behavior of the above genus expansion. 

One final word pertains to the matrix model description of strings on the resolved conifold 
via a large N duality. It was shown in [47] that there is a duality between closed and open 



topological A-model string theory on, respectively, the resolved and the deformed conifold; two 
smooth manifolds related to the same singular geometry. In the resolved conifold case the conifold 
singularity is removed by blowing up a two-sphere around the singularity; while in the deformed 
conifold case the conifold singularity is removed by growing a three-sphere around it, which 
is also a Lagrangian sub-manifold thus providing boundary conditions for open strings. As it 
turns out, the full open topological string field theory in this latter background, T*§ 3 , where we 
wrap N D-branes on the Lagrangian sub-manifold base, § 3 , reduces to SU(A) Chern-Simons 
gauge theory on S 3 |4^] , whose partition function further admits a matrix model description [|| . 
The matrix model in question, which we shall review in the next section, has a potential with 
a single minimum and no local maxima. In this paper we refer to this type of matrix models 
(which will also include the Gaussian and Penner cases) as c = 1 matrix models since, as we shall 
see, they all admit a very natural double-scaling limit to the c = 1 string at self-dual radius. 
Notice that c = 1 matrix models do not belong to the class of matrix models for which the off- 
critical instanton analysis has been carried out so far. Because understanding nonperturbative 
corrections to the topological string free energy on the resolved conifold is undissociated from 
understanding nonperturbative corrections to c = 1 matrix models, we shall consider this latter 
case more broadly in order to shed full light on this class of instanton phenomena. As such, c = 1 
matrix models is the subject we shall turn to next. 



3. c = 1 Matrix Models and Topological String Theory 

We shall now introduce three distinct matrix models, all in the universality class of the c = 1 
string, and which will be the main focus of our subsequent discussion. As mentioned in the 
previous section, one of these models is the one describing Chern-Simons gauge theory on S 3 , 
known as the Stieltjes-Wigert matrix model. Another interesting, and rather elementary, matrix 
model is the Gaussian model. Yet, we shall find that it already displays many features that will 
also appear for the resolved conifold. Finally, we also address the Penner matrix model, first 
introduced to study the orbifold Euler characteristic of the moduli space of punctured Riemann 
surfaces. These three models have been extensively studied in the literature and in the present 
section we will mostly gather some general facts necessary to obtain their topological large N 
expansions and their holomorphic effective potentials. Then, in the following section, we shall 
analyze their large N asymptotic expansions from the point of view of Borel analysis. 
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Let us begin by recalling some basic notions about matrix models (see, e.g., [49, [l8|, |], |50|). 
The hermitian N x N one-matrix model partition function is 



vol (U(N)) 



dMe-^^, (3.1) 



with vol(U(./V)) the usual volume factor of the gauge group. In the eigenvalue diagonal gauge 
this becomes 



Z = ^/n(^)A 2 We-i EC-™ (3.2, 



1=1 

where A (A) = Y\i<j (^ ~~ ^j) 1S the Vandermonde determinant. The free energy of the matrix 
model is then defined as usual F = log Z and, in the large N limit, it has a perturbative genus 
expansion 

+oo 

F = Y J 9 2 r 2 F 9 {t), (3.3) 

3=0 

with t = Ng s the 't Hooft coupling. Multi-trace correlation functions in the matrix model may 
be obtained from their generating functions, the connected correlation functions defined by 

/ 1 1 \ +0 ° 
W h (z u ...,z h ) = (Tr • • • Tr ) = V £ a+h - 2 W g , h (zi, ...,z h ; t). (3.4) 



Z\ — M z h -M 



(c) 3=0 



In particular, the generator of single-trace correlation functions is W%(z) = Nuj(z) where oj(z) is 
the resolvent, i.e., the Hilbert transform of the eigenvalue density p{\) characterizing the saddle- 
point associated to the matrix model large N limit. In the most general case, this saddle-point 
is such that p{\) has support C, with C a multi-cut region given by an union of s intervals 
Cj. At large N, the eigenvalues condense on these intervals Ci in the complex plane and one 
may interpret them geometrically as branch cuts of a spectral curve which, in the hermitian 
one-matrix model, would be a hyperelliptic Riemann surface corresponding to a double-sheet 
covering of the complex plane C, with the two sheets sewed together by the cuts Ci. 

The spectral curve, to be denoted by y(z), may be written in terms of the genus zero resolvent 
which, for a generic one-cut solution with C = [a, b], is given by the ansdtz 



1 f dw V'(w) (z-a)(z-b) 

^o{z) = — i>— W 7 V( jT, 3.5 

2t J c 27n z — w V (w — a){w — b) 

where one still has to impose that ujq{z) ~ - as z — > +oo, in order to fix the position of the cut 
endpoints 8 . The spectral curve is then defined as 



y(z) = V'(z) - 2tuj (z) = M(z) (z - a)(z - b). (3.6) 

For future reference, it is also useful to define the holomorphic effective potential, defined as 
the line integral of the one-form y(z) dz along the spectral curve, 

^h;efr(A) = f X dzy(z), (3.7) 



3 This boundary condition states that the eigenvalue density is normalized to one in the cut, J c d\p(X) = 1. 
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which appears at leading order in the large TV expansion of the matrix integral as 



N / N \ 

Z ~ / n dA, exp Yl ^h;eff (\) + •••]. (3.8) 

i=l V 5s i=l / 

Because the real part of the spectral curve relates to the force exerted on a given eigenvalue, 
it turns out that the effective potential V e g(z) = Me Vh- e s(z) is constant inside the cut C, i.e., 
inside the cut the eigenvalues are free. The imaginary part of the spectral curve, on the other 
hand, relates to the eigenvalue density as Imy(z) = 2irtp(z), thus implying that the imaginary 
part of Vh-es( z ) is zero outside the cut and monotonic inside. These two conditions guarantee 
that the eigenvalue density is real with support on C. Furthermore, as should be clear from 
the expression above, if the matrix integral Z is to be convergent a careful choice of integration 
contour for the eigenvalues has to be made based also on the properties of the holomorphic 
effective potential |2], [|. In particular, this contour may be analytically continued to any contour 
which includes the cut C and does not cross any region where V e g(z) = Re Vh;eff (z) < 0, thus 
guaranteeing global stability of the saddle-point configuration and convergence of the matrix 
integral (as MeVh ; eff(A) — > +00 at the endpoints of the integration contour). These properties 
of Vh-,efi( z ) ensure that, in the large N limit, the matrix integral can be evaluated with the 
steepest-descendent method [Q, E|. 

There are many ways to solve matrix models. In particular, |HJ] proposed a recursive method 
for computing the connected correlation functions (3.4) and the genus-g free energies, F g (t), en- 



tirely in terms of the spectral curve. This recursive method, sometimes denoted by the topological 
recursion, appears to be extremely general and applies beyond the context of matrix models; see 



[52] for a review. For our purposes of computing the genus expansion of the free energy one 
of the most efficient and simple methods is that of orthogonal polynomials [p|, which we now 
briefly introduce. If one regards 

dp(z) = e-^ V{z) ^ (3.9) 

as a positive-definite measure in R, it is immediate to introduce orthogonal polynomials, {p n (z)}, 
with respect to this measure as 

dp(z) p n (z)p m (z) = h n 5 nm , n > 0, (3.10) 

where one further normalizes p n (z) such that p n (z) = z n + • • • . Further noticing that the 
Vandermonde determinant is A(A) = detp,-_i(Aj), the one-matrix model partition function may 
be computed as 

N-l N 



Z N = n K = h$ n r»~ n , (3.11) 



n=0 n=l 



where we have defined r n = , hn for n > 1. These coefficients also appear in the recursion 
relations of the orthogonal polynomials, 

Pn+l{z) = (z + S n )p n (z) - r n p n ^ 1 (z). (3.12) 

In the large N limit the recursion coefficients approach a continuous function r n — > R(x), with 
x = Q G [0, 1], and one may proceed to compute the genus expansion of logZ by making use of 
the Euler-MacLaurin formula; see 113, ISOl for details. 
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3.1 The Gaussian Matrix Model 

Let us first focus on the Gaussian matrix model, defined by the potential Vq{z) = \z 2 . This case 
is rather simple as the matrix integral can be straightforwardly evaluated via gaussian integration, 
and the volume of the compact unitary group follows by a theorem of Macdonald [53] as 

TOl(u ' w)) = W(iv + i) ' (3 ' 13> 

where G 2 (z) is the Barnes function, G 2 (z + 1) = T{z)G 2 {z). The Gaussian partition function 
thus reads 

N 2 

Z G = ^^G 2 {N + l). (3.14) 
(2vr)^ 

The same result can be obtained with orthogonal polynomials. With respect to the Gaussian 
measure dfi(x) = e~ x dx one finds Hermite polynomials, H n (x), and for the Gaussian matrix 
model it follows 

Pn (z) = (^Y H n (-?=), h^=gy\M, (3.15) 



fWsJ V27T 

indeed reproducing the expected result for the partition function as 

N-i N-l 

z = n h n = n n! = zg > ( 3 - ig ) 

n=0 i^ 71 ) 2 n=0 

where we have also used that G 2 {N + 1) = 1^=0* n " The asymptotic genus expansion of the 
Gaussian free energy Fq = log Zq simply follows from the asymptotic expansion of the logarithm 
of the Barnes function and one obtains 

F G (i) = Wlogi-0, (3.17) 

F?(t) = -^logi + C'(-l), (3-18) 
F?(t) = = ,5 29 * g>2, (3.19) 



where £(z) is the Riemann zeta function. One immediately notices that all free energies with g > 1 
diverge when t — > 0. It is then quite obvious to consider the double-scaling limit, approaching 
the critical point t c = 0, as 



t-t c 



t -»• 0, g s -> 0, a = fixed, (3.20) 

9s 

in order to obtain the c = 1 string at self-dual radius behavior 

Finally, it is very simple to compute the one-form on the spectral curve of the Gaussian 
model 

y(z) dz = \/ z 2 - 4t dz, (3.22) 



- 18 - 



Figure 2: The Gaussian algebraic curve for values of t = —0.5, —0.1, 0, +0.1, +0.5, from left to right, 
respectively. Notice that the algebraic curve is singular for t = 0. 



3 r 




-3 |- 
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Figure 3: The real part of the Gaussian holomorphic effective potential in the complex z-plane, for 
t = 0.2. In blue is the region M.eV® eS (z) < 0, in orange the region ReV^ cS (z) > 0, and the black lines 
correspond to the Stokes lines M.eV^ efl (z) = (which also include the cut of the spectral curve). The 
white cut corresponds to the logarithmic branch cut. 

as well as the holomorphic effective potential 



Vg ee (z) = L -z^At-2t\o g j , (3-23) 

where we have normalized the result such that = 2y/i) = 0. In Figures § and || we plot 

the Gaussian algebraic curve for different values of t, as well as the real value of the holomorphic 
effective potential in the complex plane. We notice that, with an appropriate identification of 
parameters, the Gaussian holomorphic effective potential coincides with the action associated to 
the semiclassical Schwinger effect, ( p.25| ). In the following we shall further comment about this 
interesting coincidence. 
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3.2 The Penner Matrix Model 

The second example we wish to address is the Penner matrix model |34|]. First introduced to 
study the orbifold Euler characteristic of the moduli space of Riemann surfaces at genus g, with 
n punctures, it turns out that in the double-scaling limit this model is actually related to the 
usual c = 1 noncritical string theory, its free energy being a Legendre transform of the free energy 
of the c = 1 string compactified at self-dual radius J55| , |56| ] . The Penner matrix model is defined 
by the potential Vp(z) = z — logz and one may simply compute its partition function again using 
orthogonal polynomials. Indeed, one may write the Penner measure as 

dfi(z) = zss e 9s — , (3.24) 

which is, up to normalization, the measure for the generalized, or associated, Laguerre polyno- 
mials L^\x) = A e x x _a ^r (e~ x x n+a ). It thus follows for the Penner matrix model 

p n (z) = (-ir<?>! Z#/*) (£) , hF n = i- gT n\ V L + 1 + lV (3.25) 
This immediately leads to the calculation of the partition function in this model as 

A N+ id g 2 (n + i)g 2 (n + ± + i) 

*P = II hi = ^ ,} /' \ (3.26) 

n=0 ^> G 2 U + V 



where we made use of the Barnes function, satisfying 

N-l 

T(n + a + l) = ^ 

^2 



TT , x G 2 (N + (X + 1) 

Yl r n + a + 1 = 2 ^ ; — J - 3.27 



n=0 



The normalized Penner free energy .Fp = Fp — Fq = log -^7 is given by 

J-p = In 2 log <? s + - log g s - log 2^ + log G 2 ( N + - + 1 ) - log G 2 ( - + 1 ) , (3.28) 
2 5s 2 V 5 S J \g s ) 

and it admits the following genus expansion, obtained from the asymptotic expansion of the 
logarithm of the Barnes functions, 

4" (*) = \ (t + I) 2 (log (* + 1) " 5) + I (3-29) 
= -^log(t + l), (3.30) 

^ ) =2^2)(^ + 1 ) 2 " 29 - 1 )' ^ 

One immediately notices that all free energies with g > 1 diverge when t — > — 1. It is then quite 
obvious to consider the double-scaling limit, approaching the critical point t c = —1, as 

t -»• -1, 2s -> 0, /x = ^— ^ fixed, (3.32) 
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Figure 4: The Penner algebraic curve for values of t = —0.2, —0.01, 0, +0.01, +0.2, from left to right, 
respectively. Notice that the algebraic curve is singular for t = 0. 



in order to obtain the c = 1 string at self-dual radius [55] 



2-2g 



9>2. 



(3.33) 



2g(2g - 2) 

Next, let us address the large iV expansion of the Penner matrix model by making use 
of saddle-point techniques [57, 58]. This time around, the ansdtz for the large N, genus zero 
resolvent is [57] 

(3.34) 



zt \ zy ab ) 



so that its large z asymptotics, ujq(z) 
of the cut C = [a, b] to be 



+ • • • as z — > 00, immediately determine the endpoints 



1 + 2t - 2y/t(t + 1), 



b = 1 + 2t + 2y/t(t + l). 
It is now simple to obtain the one-form on the spectral curve of the Penner model 

y(z) dz -- 



2(2t + l)z + ldz 



(3.35) 
(3.36) 



(3.37) 



as well as the holomorphic effective potential 



V£ cS {z) = ^z 2 -2(2t + l)z + l + log z - log 1 - (2t + l)z+ V 'z 2 - 2 (2t + 1) z + 1 



- (2t + 1) log [z - (2t + 1) + ^-2(2^ + 1)^ + 1] + (t + 1) log (t(t + 1)) + 
+ (t + l)log4 + ivr, (3.38) 

where we have normalized the result such that Vtrcff(^) = ^- Fig ures @ an d |5| we plot the 
Penner algebraic curve for different values of t, as well as the real value of the holomorphic 
effective potential in the complex plane. The structure of Stokes lines for this potential is now 
more complicated (see, e.g., |57], |58| ) than in the familiar polynomial cases (see, e.g., P, |||). 

3.3 The Chern Simons Matrix Model 

We now turn to the Chern-Simons, or Stieltjes-Wigert, matrix model. As we previously stated 
this model is particularly interesting for its relation, via a large TV duality, to topological string 
theory on the resolved conifold [47|. The SU(iV) Chern-Simons gauge theory on a generic three- 
manifold has been realized as a matrix model in ||; see [59] for a review.. Here, we shall focus 
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Figure 5: The real part of the Penner holomorphic effective potential in the complex z-plane, for t = 0.1. 
In blue is the region ReV^. cS (z) < 0, in yellow the region ReV^. ee (z) > 0, and the black lines correspond 
to the Stokes lines ReV^ p off (z) = (which also include the cut of the spectral curve). The white cuts 
corresponds to the logarithmic branch cuts. 



on the resolved conifold case, where the partition function of SU(iV) Chern-Simons gauge theory 
on § 3 is, up to a factor, given by the Stieltjes-Wigert matrix model |j(J defined by the potential 
Vsw( z ) = \ (log^) 2 - To be precise, the Chern-Simons partition function relates to the Stieltjes- 
Wigert partition function by the simple expression Zcs = exp (— yjj (7N 2 — l)) Zsw 5 so that the 
corresponding free energies equate as 

Fcs = 4S + T2 + Fs - < 339 > 

For a review of the main features of this matrix model, including saddle-point methods and 
orthogonal polynomial analysis, we refer the reader to, e.g., [50]. 

Let us start by computing the partition function Zqs using orthogonal polynomials — as we 
shall see one may regard the Stieltjes-Wigert matrix model as a q-deformation, in the quantum 
group sense, of the Gaussian matrix model. The logarithmic measure dfi(z) = e 2sI^ losz ^ ^£ is 
well-known in the literature precisely because it leads to so-called Stieltjes-Wigert orthogonal 
polynomials, 



Pn(z) 



11 



E 

fc=0 



k(k — n) 



~2Z 



sw 



l„( n+ l) + l 




where we have introduced 



With this information at hand, one may now explicitly compute the Stieltjes-Wigert partition 
function from definition 



N-l jv N-i 

—\ 2 ^(7^ 2 -D TT r«l l 



^sw = II ^ W = f 2 ^^-^ II N 9 !- (3-42) 



,2tt> 

n=0 n=0 



A simple glance at ( [3.16|) immediately shows that, up to normalization, one may indeed regard 
the Stieltjes-Wigert matrix model as a q-deformation of the Gaussian matrix model, at least at 
the level of the partition functions. One may further define the q-deformed, or quantum Barnes 
function as 



N-l 



G q (N + l)= Yl [n] q l, (3.43) 



n=0 



so that the Stieltjes-Wigert partition function is simply Z^w = (f|) 2 q J ^^ 7N ^ 1 ^ ) G q (N + 1). 
These expressions may then be used to address the large N topological expansion of the Stieltjes- 
Wigert matrix model. Standard use of orthogonal polynomial techniques |ffq], as described, e.g., 



in []50|], yield 



9 



^o CS W = ^-^-Li3(e- i )+C(3), (3.44) 

F? S (t) = -^ + ^Li 1 (e-*)+C / (-l), (3-45) 

F cs {t) = B 2g B2 2 B 2g t . 

9 U 2<?(2 5 -2)(2 5 -2)! 2g (2g - 2)! 2g { 1 y ~ V ; 

where Li p (z) is the polylogarithm of index p, 

+oo „ 
Z 



n=l 



At genus g > 2 of course the topological expansions of Chern-Simons and Stieltjes-Wigert 



perturbative free energies coincide, F^ s (t) = i 7 "^ w (i). Furthermore, after analytical continuation 



g s — > ig s , the free energies F~ (t) coincide with the free energies of topological strings on the 



resolved conifold ( 2.42; ), once one identifies 't Hooft coupling and Kahler parameter. Finally, 



notice that all free energies with g > 1 diverge when t — >• which corresponds to e~* — > 1; with 
this second variable the natural one as the divergences are associated to the singular point of the 
(negative index) polylogarithm, Li_ p (l). It is then quite natural to consider the double-scaling 
limit, approaching the critical point e~ tc = 1, as 

e - * 1, g s -> 0, n= 6 C - 6 fixed, (3.48) 

9s 

in which case one again obtains the c = 1 string at self-dual radius 

9 2 r 2 if S (t) -> 2 J g 29 _ 2) M 2 " 2 ', 5 > 2. (3.49) 



9 Notice that, unlike in the previous example of the Penner model, in here we have not normalized the Chern- 
Simons free energy by the Gaussian free energy. 
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Finally we address the spectral curve and holomorphic effective potential for the Stieltjes- 
Wigert matrix model by making use of saddle-point techniques. With potential Vsw[z) = 
^(logz) 2 and Vg W (z) = |logz one must be a bit careful in applying (|3.5[ ) to compute the 
resolvent: indeed, the deformation of the contour around the cut, C = [a,b], must now be done 
differently due to the logarithmic branch-cut. Instead of capturing the pole at z and the pole at 
oo, this time around one captures the pole at z and the branch cut along the negative real axis 
(zero included); we refer the reader to [5(| for further details. The endpoints of the cut are 



a, b = 2e 2i - e 4 ± 2e t ^e* - 1, (3.50) 
while the one-form on the spectral curve reads 



2 1 + e~*z + J (1 + e-tz) - Az 

y(z)dz = -log v — -= dz, (3.51) 

z 2^/ z 

which coincides with the one-form log Y(Z) 4^ on the mirror curve 7i(Z, Y) = of the resolved 
conifold, written in terms of the C* variables Z = e z and Y = e y . One further computes the 
holomorphic effective potential as 

V$s(z) = -ilog 2 z + log 2 C-2(logeiog(l-e- t O+Li 2 (l-0 + Li 2 (e-*e)) -14) 
= -~ log 2 z + log 2 £ - 2 ( log e log(l - e-*f ) + Li 2 (e"*e) - Li 2 (£) - 

-Iog(l-01oge + y)-^o, (3-52) 
where equality holds due to the Euler's reflection formula for dilogarithms 

Li 2 (0 + Li 2 (l - = - log (1 - log £■ (3.53) 
b 

In here we have set 



1 + e~*z + J(l + e _t z) - 4z 
= — 2 (3-54) 

to simplify notation 10 and we have defined 

V = -- log 2 (2 e 2 ' - e* + 2 e'f yV - l) + log 2 (V + el \/e* - lj - 

-2 log (V + e5 yV - l) log (-e~^vV - l) - 2 Li 2 (l - e f - e5 yV - l) - 
-2Li 2 (l + e"Ve* - l) , (3.55) 



to ensure that the result is normalized such that ^vJr(^) = ^ n Fig ures |6| and ^ we plot the 
Stieltjes-Wigert algebraic curve for different values of t, as well as the real value of the holo- 
morphic effective potential in the complex plane. Given that Li 2 (z) is the standard dilogarithm 
function, with its intricate branch structure, it is not too hard to realize that the structure of 
Stokes lines of the present effective potential is now much more complicated than usual. 



Using this variable, the spectral curve is also compactly re-written as y{z) = -j log 
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c 



c 



Figure 6: The Stieltjes-Wigert algebraic curve for values of t = —0.1, —0.01, 0, 
right, respectively. Notice that the algebraic curve is singular for t = 0. 



-0.01, +0.05, from left to 




Figure 7 : The real part of the Stieltjes-Wigert holomorphic effective potential in the complex z-plane, 
for t = 0.1. The black lines correspond to the Stokes lines ReV^ s (z) = (which also include the cut of 
the spectral curve). The white cuts and regions correspond to the logarithmic and dilogarithmic branch 
cuts. Because of the choice of principal sheets in Mathematica the colored regions are now not so clear. 
Akin to the Penner model, the region in yellow, to the right of the vertical black lines, and the region in 
blue, inside the "closed bubble", have feV^^z) > 0, in the principal sheet. The rest is ReVSf(^) < 0. 



3.4 Double Scaling Limit and c = 1 Behavior 

We have just seen that the Gaussian, Penner and Chern-Simons free energies admit rather simple 
double-scaling limits to the c = 1 string at self-dual radius. In the Chern-Simons case this relates 



to our earlier discussion in section |2.3| , where we pointed out that, at the conifold point of moduli 
space, the A-model may still be characterized by its critical behavior in the double-scaling limit 
(2.32). Free energies of the topological string reduce, in this situation, to free energies of the 
c = 1 string. Let us now briefly discuss, in the example of the Chern-Simons matrix model, how 
one may also study the destiny of the open string sector, i.e., of the matrix model correlators 
Wh{pi, ■ ■ ■ ,Ph) introduced in ( |3.4[ ), in this c = 1 double-scaling limit. 
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Introducing a parameter £ as 

e -t = 1 - C, (3.56) 

the conifold point of the Chern-Simons model is thus located at £ = 0. The expansion of the 
branch points of the spectral curve ( |3.51 ) near the conifold point yields 

a,b = 1±2C3 + ••• , (3.57) 

in which case it is natural to scale also the z variable in the spectral curve as 

z = 1 + £5 s (3.58) 



in order to appropriately zoom into the critical region. In here, s is the double-scaled open 
coordinate. In these variables the Chern-Simons one-form y(z) dz, at criticality, scales to 



y(z) dz -> y(s) ds = ( \As 2 - 4ds (3.59) 

which one immediately recognizes as the one-form of the Gaussian matrix model. The interesting 
point is that, in the same variables, also the two-point correlator Wo(p, q) reduces to the Gaussian 
one 

Wo 2 (p,q)dpdq^ - . 1 N9 ( - St ~ 4 -l)dsdt. (3.60) 



In fact, there is a property of the topological recursion, proved in |5T|], which states that one may 
either first compute matrix model amplitudes and then take their double-scaling limits, or else 
recursively compute amplitudes directly from the double-scaled curve, the result being the same 
(i.e., the operations commute). As such, in the c = 1 double-scaling limit Chern-Simons open 
correlators will all reduce to Gaussian open correlators 

W$ (*!,..., z h ) dz x --- dz h -> ( 2 - 2 °- h Wg, (si, . . • , 8 h ) dsi • • ■ da fc . (3.61) 

Now recall that open topological string amplitudes may be computed from the matrix model 
correlators W 9t h(zi, . . . ,Zh) as |n], 62, 10, |Tl| 

i \ f Pl f Ph 

A ( h 9 \ Pl ,...,p h ) = J ■J dz 1 ---dz h W g>h {z 1 ,...,z h ), (3.62) 

where the {pi} are the open string parameters which parametrize the moduli space of the brane. 
As such, ( |3.61| ) shows how, near the conifold point, open amplitudes of topological strings on the 
resolved conifold reduce to Gaussian amplitudes. This is actually generic for topological string 
theory near the conifold point [12). We shall now relate these Gaussian amplitudes with open 
amplitudes in the dual c = 1 model. 

We first need to recall some results in non-critical strings, holographically duals to matrix 
models. We are interested in minimal models obtained by coupling 2d gravity to minimal (p, q) 
matter models, with central charge c p>q = 1 — 6(p — q) 2 /pq- The coupling to gravity leads to the 
appearance of the Liouville field, eft, with world-sheet action 

s L = f^-Va (V + QR0 + 4 W 2 ^) , (3-63) 

where /iL is the bulk cosmological constant. The central charge of the Liouville sector is cl = 
1 + 6Q 2 , and the parameter b above relates to the background charge Q as Q = 6+1/6. The 



- 26 - 



bosonic string requirement that the total central charge of Liouville theory plus minimal matter 
equals c = 26 eventually fixes b = ^/|. There are two distinct types of boundary conditions 

in Liouville theory ]63| , |64[ . There is a one-parameter family of Neumann boundary conditions, 
the so-called FZZT branes, parameterized by the boundary cosmological constant /^b, usually 
expressed in terms of a parameter s as 



A*B = A / . ^ L , 9 , cosh (nbs) . (3.64) 
y sin (7ro z j 

Besides FZZT branes, there are also ZZ branes, associated to Dirichlet boundary conditions. 
These correspond to a two-parameter family, parameterized by the pair of integers (m,n), and 
are localized at <p = oo. At the quantum level FZZT and ZZ boundary conditions, or, respectively, 
the (B s \ and (J3( mn \| boundary states, are related as |65|, ^6| |67fl 

i B (m,n) I = { B s(m,n) I ~ ( B s(m-n) I , with s(m, n) = i + bnj . (3.65) 

Both types of branes have been given a geometrical interpretation in terms of a complex 
curve, in |37|]. This is accomplished by introducing the variables 

® -7FZZT . , fts 



x = u B ~ cosh (vr6s) , y = - — Z bLLl ~ sinh — . (3.66) 

OfXB \ J 

Considered as complex variables, the coordinates {x, y} define an algebraic curve F(x, y) = 
embedded into C 2 , which is identified with the spectral curve of the dual matrix model {i.e., a 
double-scaled hermitian one-matrix model) j67j. The FZZT brane disk partition function may 
be equivalently written as the line integral of the one-form y Ax as 

Z FZZT M= f^dxy. (3.67) 



Analogously the /i-point matrix model correlators (3.4) are identified with open amplitudes with 
FZZT boundary conditions, e.g., the two-point function Wo,2(p,<l) is identified as the annulus 
amplitude for FZZT branes and so on. The ZZ brane disk partition function is instead defined 
as the line integral of y dx over a closed contour 

Z (m,n)=<f dxy, (3.68) 

J 7m,n 

where 7 m) „ is a non-contractible contour conjugate to a "pinched cycle", starting and ending at 
the singular point X( m , n ) = x (s(m,n)) and y (m) „) = y (s(m,n)) |37|. 

The case of c = 1 is a bit more subtle since one has to consider the singular b — >• 1 limit. It 
is first necessary to introduce the renormalized couplings 

Hc=l = lim (tt (l - b 2 ) hl) , /Ub,c=i = lim (it (l - b 2 ) /is) • (3.69) 

Furthermore, an appropriate subtraction is required in order to define the FZZT disk partition 
function. This may be expressed in terms of the one-form w(s) d//B,c=i(s) with 

/a vFZZT 4 \ 

w(s) = lim ^ + - Z D /is , (3.70) 
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where Zd is the disk partition function in the c = 1 CFT. The relevant c = 1 curve then reads 



x(s) = A*B,c=i(s) = VA*c=i cosn (ns) , y(s) = u;(s) = -D^/j, c=1 its sinh (tts) , (3.71) 



where D is some constant. The identification of the curve ( |3.71| ), arising from CFT considerations, 
with the curve of the dual matrix model is another delicate point. Here, the relevant matrix 
model is a double-scaled version of Matrix Quantum Mechanics (MQM) with a Sine-Liouville 



perturbation. It is know for quite some time [49] that the singlet sector of this MQM can be 
reduced to a system of free fermions, in an inverted harmonic oscillator. In the semiclassical limit 
the ground state of this system is completely determined by the shape of the Fermi sea, which 



can be parameterized in terms of an uniformization parameter r as [40] 



x(t) = y2/i cosh(r), y(r) 



2fi sinh(r), 



(3.72) 



where p denotes the Fermi level. In analogy with the c < 1 case, one would like to identify the 
above MQM curve with the CFT curve ( p. 71 ). However, these two curves are clearly distinct. A 

where it was proposed that one should 



70] 



solution to this puzzled has been offered in [68, 
instead identify w(s) with the resolvent, rather than directly with the spectral curve of the dual 
matrix model. The spectral curve can then be extracted, following a very standard matrix model 
procedure, from the discontinuity of w(s), 



P(s) 



1 
2vri 



(w(s + ie) — w(s — ie)) = —D^J p c= \ sinh (tts) . 



Clearly the new CFT curve, defined as 



>( a ) = Mb,c=i(s), y(s)=p(s), 



(3.73) 



(3.74) 



agrees with the one-matrix model spectral curve ( 3.72j ) after an appropriate identification of 
parameters. The FZZT brane partition function can then be obtained from the line integral 
of the one-form y(r) d r x(r) dr, while the ZZ brane partition function can be defined by the 



following closed integral on the MQM curve 6S] 



Z (n,l) 



dx y 



In 



dr d T x(r) y(r) = 27rn/i, n G Z, 



(3.75) 



corresponding to a (n, 1) ZZ brane partition function. Indeed, it has been shown that only an 
one-parameter set of the c = 1 ZZ branes may be identified in the dual MQM J71| |. 

Let us further notice that the above matrix quantum mechanics spectral curve, ( gTrD , with 
the uniformization parameter r, is just an infinite covering of the hyperboloid 



y 2 = 2m, 



(3.76) 



which is precisely the spectral curve of the Gaussian matrix model. In particular, this explains 
how open matrix model correlators of the Gaussian model, W^ h , get identified with D-brane 
amplitudes with FZZT boundary conditions in the c = 1 model at self-dual radius. Indeed, in 
[72] it was checked that the double-scaled Gaussian correlators are related to macroscopic loop 
operators in the c = 1 theory. Finally, the limit ( p. 61 ) shows that topological string amplitudes 
with toric-brane boundary conditions reduce to c = 1 amplitudes for FZZT branes. Hence, and 
as already pointed out in a related context in [jn|, toric branes reduce to FZZT branes in the 
double-scaling limit, at the conifold point. 
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4. Nonperturbative Effects, Large Order and the Borel Transform 

We may now turn to the study of the asymptotic perturbative expansions for the free energies of 
the matrix models and topological strings we are interested in. In particular, we shall perform 
a detailed Borel analysis of each case, and thus understand what type of nonperturbative effects 
control the large-order behavior of the distinct perturbative expansions. 



4.1 The Gaussian Matrix Model and c = 1 Strings 

Let us begin with the Gaussian matrix model. The genus expansion of its free energy, (3.19), is 



clearly an asymptotic expansion withF 5 G ~ (25-3)!, given the growth of Bernoulli numbers as 
B2g ~ (2<?)!- Recalling our discussion in section g, we may then consider the Borel transform of 
the divergent Bernoulli sum — i.e., restricting to genus g > 2 — and obtain 



This function has no poles on the positive real axis, for real argument (the genus expansion (3.19) 
is an alternating series). As such, one can define its inverse Borel transform 11 

M9s) = A f ^ - ( ■ V,,. x " I " ) - + ^ ! < ( 1-2) 




4 Jo s ysinh 2 (ffs) 

providing a nonperturbative completion for the asymptotic expansion of the free energy in the 
Gaussian matrix model. It is quite interesting to notice that, upon the trivial change of variables 
s — )■ a = || s, this expression precisely coincides with the one-loop effective Lagrangian for a 
charged scalar particle in a constant self-dual electromagnetic field (of magnetic type) introduced 
in section |2.2j . Comparing with ( p. 20 ) we see that in here 



2eJ7 1 (A 

m z iV 



If one instead considers imaginary string coupling, g s = ig s , the asymptotic expansion ( [4.2|) 
will coincide with the one-loop effective Lagrangian corresponding to a self-dual background of 
electric type, which is exactly the same as that for c = 1 strings at self-dual radius. This time 
around the perturbative series is not alternating in sign, and the Borel integral representation 



1 f + °° da ( 1 1 1\ _2ta 

F c=1 (g s ) = - / — — a-r e - 4.4 

4 Jo o" Vsm^cr a z 3/ 



has an integrand with poles on the positive real axis, in principle leading to ambiguities in the 
reconstruction of the function, as discussed in an earlier section. However, we may now use the 



analogy of this expression to the results in section 2.2 in order to use the unitarity prescription to 
perform an unambiguous calculation — which basically yields an ie prescription which reduces the 
imaginary part of the integral to a sum over the residues of its integrand. The nonperturbative 



11 Notice that since ~ (2g — 3)! the inverse of the Borel transform will now have an extra factor of ^ with 
respect to the definition of section H, which dealt with asymptotic growths of the type ~ (/3n) ! . 
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imaginary contribution to the above free energy is thus simple to compute as 

ImF c= i {g s ) = - > f — - — -2 j ~ q e Ss 

4 „ J titt 2vri cr V sir a a z 3 ' 



[2 



n=0 Jn7r 



, g (? + |W ^. ( , 5) 



As expected from the discussion above, this formula precisely matches with the Schwinger result 
in a self-dual background, expressed in ( p. 22 ). It can also be obtained from the "alternating" 
result ( [4.2[) by analytic continuation and contour rotation. Furthermore, as discussed in section 



2.1, it follows that the discontinuity of the free energy across its branch cut consists of an instanton 
expansion given by 

Disc iUC*) = £ ( V + I) ^ = + ^Gfc) + • • • • ( 4 - 6 ) 

We may now relate this instanton series to the full c = 1 perturbative expansion, ( 2.33 ), 
by means of the Cauchy formula fl2.8|) . One first observes that the integral over the contour 
at infinity in ( |2.8[ ) has, in here, no contribution, since the Barnes function is regular at infinity 
(see, e.g., JT3|). As such, the dispersion relation fl2.£|) reads 13 , after power series expansion of the 
integrand's denominator, 

v^\-^ qt f dz I 2irtn\ _-^tn 



- gists <r™ - g «j#<«-<* -> (f )'".«. 



where we used the definition of the Riemann zeta function as 

+00 



1 n 

n=l 

Notice that, from the first to the second line, we truncated k = from the A; sum. Indeed, 
for this particular value of k the integral would require regularization. However, this would only 
contribute to terms at genus zero and one, which we are not considering here in any case. As such 
we shall simply truncate the k = contribution from the sum, without the need to regularize 
the divergence, and focus on the genus g > 2 contributions. In this way, if in the above formula 
for F c= i(g s ) we further relate the Riemann zeta function to the Bernoulli numbers via 



C(2n) = (-l)-(" +1 )^^B 2n , (4.9) 



it immediately follows 



-00 



\B 2 ,\ fS,\ 2 '-' 2 



9=2 



12 Notice that the pole at a = has vanishing residue. 



Recall that for matrix models and strings one uses z = g a ; see section 2.1 
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which is indeed the c = 1 perturbative expansion, for genus g > 2, with fj, = 



In some sense, ( 4,10| ) takes us back to where we started the discussion, i.e., the Gaussian 



matrix model perturbative series. Indeed we have seen that the alternating Gaussian perturbative 
series admits a simple Borel transform, which may be inverted unambiguously to provide a 
nonperturbative completion of the theory. Upon "Wick rotation" of the coupling constant, this 
completion also describes the non-alternating c = 1 string theory alongside with its instanton 
effects (obtained in a fashion very similar to our earlier discussion of the Schwinger effect). Of 
course that a key aspect of this analysis is the fact that the integral representation of the free 
energy, provided by the inverse Borel transform, precisely coincides with the nonperturbative 



integral formulation of the c = 1 theory put forward in [25, 43|. One may thus consistently pick 
either starting point and obtain the very same results. 

To end our analysis, we shall now address the large-order behavior of perturbation theory 
and see that it is controlled — as expected — by one-instanton contributions, i.e., by the closest 
pole to the origin in the complex Borel plane. If one considers the first term in the instanton 



expansion (|4.6|) and, following the discussion in section 2.1, one sets 



F^(9s) = -^(t + ^)e-^, (4.11) 
g s v 27T/ 



a comparison with ( 2.1 0| ) immediately yields 



A = 2nt, b = -l, Fi 1] =t, F 2 (1) = -^, (4.12) 

thus identifying the instanton action, the characteristic exponent, and the loop expansion around 
the one-instanton configuration. In fact, it is rather interesting to observe that in this situation 
the loop expansion around each ^-instanton sector is finite. This is quite unusual; typically the 
^-instanton loop expansion is itself asymptotic, with its large-order behavior being controlled 
by the (I + l)-instanton configuration. What we observe is that in this case only the zero- 



1 



instanton sector displays non-trivial large-order behavior. Now, with the identifications ( 4.12 ) 
the large-order equation ( 2.1 1| ) implies 



} (2 g -l)(2 ff 3)! 

Checking the expected large-order behavior in this case is straightforward, as one simply needs 
to use the standard relation 

= 2 ~i$ c{29) = j^( 1+ e ~ w + e ~ 25los3 + ■ ■ ■ ) (4 - 14) 

in the asymptotic series ( |4. 10|) and the above fl4.13|) immediately follows; exponentially suppressed 
contributions in ( 4.14D are not contributing to the large order of the zero-instanton sector. If one 



is instead interested in the large-order behavior of the Gaussian matrix model, one essentially 
just needs to use the analytically continued instanton action A = 27rit and everything else follows 
in a similar fashion. 

4.2 The Penner Matrix Model 

Having worked out the Borel analysis of the free energy in the Gaussian matrix model, which 
essentially reduces to the Borel analysis of the logarithm of the Barnes function G2{z), we have 
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all we need in order to write down the nonperturbative part of the free energy of the Penner 
model (3.31). In fact the whole procedure is essentially the same as before and we shall leave 
most calculations to the reader. The Borel transform is now 
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(4.15) 



It should be simple to spot the similarities to the Gaussian case, as the free energy of the 
Penner model may be written in terms of Barnes functions as in (ET2Sj) . As such, the total 

discontinuity (for g s — > ig s ) is now given by the sum of the discontinuities of logG72 + 1 
and log G2 ( — + 1 J , which yields 
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(4.16) 



This is quite simple to obtain by following a procedure identical to what we used in the Gaussian 
case, but starting from the above Penner Borel-transform. Notice that in this case we have two 
sets of nonperturbative contributions, with instanton actions 2n (t + 1) and 2tt, respectively. The 
large-order behavior of the theory is controlled, as usual, by the closest pole to the origin in the 
Borel plane. For t > the relevant pole is located at 2ir; however, close to criticality t — > — 1, 



the first instanton tower in ( 4.16 ) is the relevant one 



4.3 The Chern Simons Matrix Model and the Resolved Conifold 

We may now turn to the Chern-Simons matrix model, holographically describing topological 
strings on the resolved conifold. The genus expansion of its free energy, ( |3.46D , is asymptotic, 
and in here we wish to analyze this divergent series from the viewpoint of Borel analysis, as we 
did earlier with both the Gaussian and the Penner matrix models. 

Let us start with the Chern-Simons genus expansion ( p. 46 ) where, for the moment, we 
drop the constant map contribution and focus on genus g > 2. This will allow us to better 
understand the divergence of the series arising from the term with both a Bernoulli number and 
a polylogarithm function contributions. One has in this case 14 



Fcs(g s 
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2g{2g-2)\ 3 " 2 ^ 6 > j^ 9s 2g (2g - 2) ^ ( t + 2nim) 
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(4-17) 

where we have used an integral representation of the polylogarithm in terms of a Hankel contour, 

(4.18) 



re-written as a sum over residues [74], to express 
Li 3 - 2g (e-*) = r (2g - 2) £ 



(t + 2^i77i) 2 ^ 2 ' 



9>2. 



Now, since the Bernoulli numbers grow as Big ~ (2g)\ and the polylogarithm functions behave, 
in worse growth scenario 15 , as lirniji-^o ^h-2g ( e_i ) ~ r (2g — 2) t 2 ~ 2g , this series is asymptotic 



'Notice that the m = contribution, in the sum in the second expression, equals Y^l 



,2g-2 



which 



^9=2 2 3 (2 9 -2) V t ) 

is of course the Gaussian free energy at genus g > 2. This is a consequence o f wo rking with a Chern-Simons free 



energy which is not normalized against the Gaussian free energy; see section 3.3 

15 At large t the polylogarithm's growth is not factorial in genus, as one has limR e t-»±oo L13- 



29 



= 0. 
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and, like in Gaussian and Penner models, its coefficients grow factorially as (2g — 3)!. In this 
case one is led to the Borel transform 



+00 TpCS( f \ +00 R t2o-2 

=£ ^ +g ±^o!_i 1 y (4 , 9) 

This function has no poles in the positive real axis, for real argument. This is expected since we 
started off with an alternating sign expansion and, in this case, we may define the free energy 
via the inverse Borel transform 

~ , . 1 ^ f + °° ds ( 1 /2(t + 2vrim)\ 2 1 l\ , . . 

Fcs(g s ) = -7 E / — 77 7 " - + o e • 4 - 20 

^Jo s ^sinh 2 ( 5( ^ S ) V 9s J s 2 3) 

This provides an unambiguous nonperturbative completion for the asymptotic expansion of the 
free energy in the Chern-Simons model. Interestingly enough, and like in previous examples, for 
each distinct m a trivial change of variables turns the corresponding expression into the one-loop 
effective Lagrangian for a charged scalar particle in a constant self-dual electromagnetic field (of 
magnetic type) introduced in section |2.2| , the sum over all integer m thus corresponding to a 
sum over an infinite number of Lagrangians of this type. If we further report to the discussion 



in section 2.4, we see that our Borel resummation is essentially (up to analytic continuation, as 
we move between alternating and non-alternating perturbative series) equal to the GV integral 
representation of the free energy of topological strings on the resolved conifold. 

Let us thus consider the case of the resolved conifold in greater detail, which is obtained by 
the simple analytic continuation g s — > ig s , and corresponds to the electric version of the above 
result. In this case the free energy perturbative series is non-alternating and the Borel transform 
has poles on the positive real axis, making the reconstruction of the free energy possibly affected 
by nonperturbative ambiguities, which we may, however, understand in the computation of the 
imaginary part of the integral, 



conif 



2 s =AZ^ —( — 2 " o e SS 4 - 21 ) 



Moreover, since (4.21) agrees, after a simple change of variables, with the GV integral represen- 
tation ( 2.35 ), at least for genus g > 2, we have a physical interpretation for the nonperturbative 
terms we find: as observed earlier, in section |2,4[ , the imaginary part of the integral will compute 
the BPS pair-production rate in the presence of a constant self-dual graviphoton background. 

The imaginary part of the integral ( [4.21| ) may be computed by the use of the unitarity +ie 
prescription, yielding a sum over residues of the integrand. Equivalently it equals one-half of 
the integral over the whole real axis (the imaginary part is symmetric), which may be computed 
by closing the contour on the upper half of the complex plane, thus enclosing the poles of the 
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Figure 8: Poles in the complex Borel plane, for the case of the resolved conifold free energy. The poles 
in red coincide with the spurious c = 1 string contribution, arising from the Gaussian normalization. 



hyperbolic sine; see Figure [||. It follows 

s 7T^v-^/dCTl/l 1 1\ 2(t+27rim) 

1 f 2tt (t + 27rim) g s \ / 2-7T (i + 27rim) n 

~ 47ro s ^-J At V « n 2 / 6XP I 5s 

Mi n=l meZ v / \ us 

One observes without surprise that this formula matches (an infinite sum of) the Schwinger result 
in a self-dual background. The discontinuity of the free energy across its branch cut is thus given 
by the instanton expansion 



Making use of 
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this discontinuity may also be written as 
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(4.25) 

As we did in the previous cases, we may now use the Cauchy formula to relate this instanton 
series to the perturbative expansion of the resolved conifold's free energies. Once again, the 
integral over the contour at infinity in ( p.8| ) has no contribution (see the appendix), and the 
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dispersion relation thus reads, successively, 
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where we made use of 



/ " dx f(x) 5 (g(x)) = T^r v (4-27) 
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with Xi the real simple roots of g(x); of the definition of the polylogarithm of index p 

+00 n 

Li ^) = E^= ( 4 - 28 ) 
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and the definition of the Riemann zeta function alongside with its relation to the Bernoulli 
numbers. This result shows that the instanton expansion ( 4.23| ) indeed has enough information 



to completely rebuild the full free energy perturbative expansion for topological strings on the 
resolved conifold, at genus g > 2 (at genus and 1 one still needs to take into consideration the 
relation between Chern-Simons and Stieltjes-Wigert perturbative free energies, and additional 
regularizations may be needed as, e.g., in the Gaussian case). 

One may further show that, in particular, the closest pole to the origin in the Borel complex 
plane controls the large-order behavior of perturbation theory, corresponding to the familiar 
one-instanton contribution. A glance at Figure || makes it clear that the closest pole to the 
origin corresponds to the one-instanton contribution of the c = 1 string. This, of course, is due 
to the fact that the Chern-Simons free energy in ( p.46| ) is not normalized by the Gaussian free 
energy. This is easily corrected by considering, in the following, the normalized free energy 

-^coniffe) = F conil (g s ) - F c=1 (g s ). (4.29) 

This will guarantee that our large-order tests will precisely look at the true "resolved conifold 
contribution", without being plagued by ghost effects due to the Gaussian measure. In this 
normalized case, and as is simple to check by looking at Figure |8] again, we have two complex 
conjugate poles equally distant from the origin, at £ = 2n {t ± 27ri) /g s , and we have to consider 
the contributions from them both. Let us consider the first terms in the instanton expansion 
(|423|) 



and thus define, following section 24 
A± = 2vr (t ± 2m) => 
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A± = \A\ e ±i9A = 2tt\A 2 + 4vr 2 exp f±i arctan^ ,(4.31) 
F$ = I e ±Wp = Vt 2 + 4tt 2 exp f ±i arctan y ^ , (4.32) 

(4.33) 



Again, as in the previous examples, the loop expansion around each f-instanton sector is finite. 
Thus, also for the resolved conifold only the zero-instanton sector displays non-trivial large-order 
behavior. In particular, the large-order equation ( p. 11 ) implies 
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The check of this behavior is now harder than before, and we shall need to perform numerical 
tests to confirm its validity In this sense, one constructs the test ratio 
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(4.35) 



where equality holds at large g due to ( 4.34| ). In Figure ^ we can see that numerical analysis 
undoubtedly confirms our prediction. 

The analysis so far has focused only on the contribution to the resolved conifold free energy 
arising from the D2D0 bound states of branes. As is clear in ( 3.46| ), to this term one must still 
add the contribution from bound states of DO-branes, i.e., the contribution of constant maps. 
As it turns out, in this case the discontinuity is obtained from ( 4.23j ) by simply setting t = in 
that expression, i.e., 
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A procedure that should be familiar by now yields back the perturbative expansion ( 2.3§| ) or 
( 3.46[ ) via the Cauchy formula and an integration of the above discontinuity over the free energy 
branch cut: 
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where we have used familiar properties of the zeta function and Bernoulli numbers, including 

(4.38) 



(-D n 



n + 1 



- 36 - 




-1.0 



Figure 9: Test of large-order behavior for the resolved conifold. The plots represent the test ratio R g (t), 
in red, versus the expected behavior cos((2g — 2)arctan^), in blue, as a function of genus g, and for 
t = 20, 40, 60, 80, from left top to right bottom, respectively. The matching at high genera is evident. 



The last thing we want to show is that, as expected, it is the closest pole to the origin in the 
Borel complex plane that controls the large-order behavior of the theory. Repeating our earlier 
discussion we find the one-instanton contribution 
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In this case the large-order equation ( p. 11 ) yields 
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To check that this is the right answer all one has to do is to use 
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Notice that the above exponentially suppressed contributions do not contribute to the large order 
of the zero-instanton sector. 



5. Stokes Phenomena and Instantons from Hyperasymptotics 

Having understood the Borel analysis of topological strings and c = 1 matrix models, benefiting 
in this course of the identification of instanton effects in these models, we shall now make a 
brief detour into the realm of hyperasymptotics, as first introduced in |T5| (see, e.g., |76| for a 
review), i.e., a series of techniques to refine optimally truncated asymptotic expansions by the 
inclusion of exponentially small contributions. In particular, our focus will concern hyperasymp- 
totic approximations for integrals with saddles [ff^ , ff8| , the prototypical example for problems 
dealing with the calculation of partition functions or free energies. In this case, the exponentially 
suppressed contributions arise from saddles other than the one chosen in the steepest-descent 
asymptotic approximation. The main interest of this analysis for the present work is that the 
Stokes phenomenon — certain "discontinuities" which we shall explain below and later relate to 
instanton effects — is automatically incorporated into the hyperasymptotic scheme. 

Suppose one wants to use the method of steepest descents in order to find an asymptotic 
expansion, as \n\ — > oo with k = \k\ e lS , of the one-dimensional "partition function" 

Z(k) = [ dze- KW{z \ (5.1) 



Jc 

where C is a contour we specify below. A typical calculation goes as follows: one begins with 
the calculation of the saddle points of the "potential function" W(z), computed as the set of 
points {zk}k=i,2,- such that W'(zk) = 0; then, chosen a reference saddle-point z n , the contour 
of integration C is deformed to the infinite oriented path of steepest descent through z n , which 
we shall denote by C n (6). This contour is defined as 

Em [as (W(z) - W{z n ))] = (5.2) 

and with k(W(z) — W{z n )) increasing away from z n . This immediately implies that the phase 
of W(z) — W(z n ) must equal —9 + 2-Km, m G Z. We then introduce the "partition function" 
Z n (n) evaluated on the n-th saddle 

Z n ( K ) = -L e - KW ^Z n ( K ), (5.3) 



K 

Z n ( K ) = / dze-< w ^~ w ^\ (5.4) 

As we shall see, Z n {n) will display Stokes phenomena in the form of a discontinuity associated to 
a jump in the steepest-descent path whenever it passes through one of the other saddles, k ^ n. 
The integral fl5.4j ) can be evaluated via the steepest-descent method and one obtains a function 
of k for each saddle, n, given by a series in negative powers of k 



9=0 



with (see Jf7, [Tg] for details) 
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Figure 10: On the left, the complex k plane, showing a region — a nmi < 9 < —a nm2 . On the right, the 
equivalent region in z space, with saddle z n , adjacent saddles z mi and z m2 , and respective steepest-descent 
contours through z n (in red) hitting the adjacent saddles. In blue we plot the adjacent contours. 



The series in (^^) is asymptotic. The point of view of |77J is to understand this divergence as a 
consequence of the existence of other saddles {zk^ n }, through which C n does not pass. Because 
one is free to choose the reference saddle n at will, all possible asymptotic series are thus related 
by a requirement of mutual consistency, also known as the principle of resurgence: each divergent 
series will contain, in its late terms, and albeit in coded form due to their divergent nature, all 
the terms associated to the asymptotic series from all other saddles. 

Another important point concerning the asymptotic series ( |5.5D dwells with the fact that 
this expression only holds in a wedge of the complex ft-plane, i.e., for a restricted range of 9, a 
property which is associated to Stokes phenomena. Suppose that in the above set-up, and once 
(|5.5| ) has been computed, we start varying 9 in such a way that we always choose the contour of 
integration to be the steepest-descent through the saddle z n . As it turns out, this is a continuous 
process only for a finite range of 9: indeed, one faces a discontinuity if 9 reaches a value such 
that the contour of integration passes through a second saddle, z m . This will happen when 9 
reaches the value — a nm , with 

-o-nm = ~ arg (W(zm) - W{z n )) . (5.7) 

For this value of 9 the steepest-descent contour will change discontinuously and exponentially 
suppressed contributions to fl5.4|) will "suddenly" become of order one. In this case, and in order 
for the steepest-descent contour to go through a single saddle, one must restrict 9 to an interval 
-from < 9 < -cr nTO2 , where z mi and z m2 are saddles adjacent to z n , i.e., saddles which may 



be reached from z n through steepest-descent paths 16 . This is illustrated in Figure 1C. One of 



the goals of hyperasymptotics [75] is to deploy resurgence in order to better understand Stokes 



phenomena, and this is the aspect we shall be mostly interested in. 



16 A saddle z m is said to be adjacent to the saddle z n iff there is a path of steepest descent from z n to z m , i.e., 
z m will be adjacent to z n whenever 9 = — a nm and thus arg (W(z) — W(z n )) — a nm (naturally this is also the 
condition that defines the Stokes lines for Z n (n)). One similarly defines the adjacent contour through the adjacent 
saddle as the steepest-descent contour C m (— <J nm ), through z m . 
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Let us make these ideas more precise. In hyperasymptotics one begins with (optimal) trun- 
cation of the asymptotic series fl5,5|). In this case, 



N-l r , , 
9=0 K 

where 1Z^\k) is the remainder associated to the finite truncation. The main contribution of the 
hyperasymptotic calculation in [^] was to produce an expression for the reminder which led to 
an exact resurgence formula for ( |5.5| ), and we shall now present these results. Let us first define 
the singulant, for every adjacent saddle z m , as 

Wnm = W(z m ) - W(z n ) = \Wnm\ ^ ■ (5.9) 

In this case the remainder term (see |77| , |78[| for details) can be expressed as a sum over integrals 
through all adjacent saddles to z n , {z m }: 

ni " ){K) = ^ri ^ E / dr)^-je-^-Z m ( V ), (5.10) 

where the Z m {r\) are defined on the adjacent contours C m ( — 0"nm)' There are several interesting 
points to this formula. First, it provides an exact and explicit expression for the reminder and 
one now explicitly sees that the divergence of the asymptotic series (|5.5|) is directly related to the 



existence of adjacent saddles. Second, inserting the above expression back in (5^), one obtains 
the exact resurgence formula 

N ~ 1 r (n\ 1 1 ^ /-oo-e- i<T ™ n N - 1 

= E^ + ^E/ ^f-r ^ 2 rn(,), (5.11) 



which is the basis for the hyperasymptotic analysis of [77\: indeed, each Z m {rf) in the above 
integrands may itself be expanded as an asymptotic series leading, via iterations of the above 
formula, to exponentially improved asymptotic results for the original Z n {n) (in the sense that 
the error associated to the approximation is reduced from polynomially small to exponentially 
small 17 ). Third, as the resurgent formula holds for any N, one may write it down for N = 0, 
obtaining either the resurgent expression [ f78| ] 

Z ^ K ) = E / dr ? — 1~ Zm{rql (5 - 12) 

which leads to interesting functional relations in selected examples; or the (formal) resurgent 
relation |77j 

m h=0 vv nm 

which expresses the late terms (3 > 1) of the asymptotic series at a given saddle as a sum over 
the early terms of the corresponding asymptotic series at the adjacent saddles. In particular, the 



For a recent discussion in the field theoretic context see 
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leading contribution arises from the adjacent saddle m* with smallest singulant, i.e., to leading 
order one obtains 

(g-Dl 



C fl (n) 



W 9 , 

' " nm* 



CoK), 



(5.14) 



which makes manifest the characteristic factorial behavior of asymptotic series. Fourth, and 



finally, the resurgence formula ( 5.11 ) precisely incorporates Stokes phenomenon 77]; the appear- 
ance of suppressed exponential terms as the steepest-descent contour C n (9) sweeps through one 
of the adjacent saddles, m (see Figure |i"T| ), This, as we mentioned, will happen as 9 crosses the 
Stokes line C n (—a nm ), in which case the asymptotic expansion (|5.5[) will have a discontinuity 



(5.15) 



AZ n { K ) = Z n (\k\ e i (-' Tnm+0+ )) - Z n (\k\ e i (- anm+0 )) / 0. 



It is not too hard to compute the precise value of this discontinuity straight from the resurgence 
formula for the remainder, fl5.ll ). One obtains, without surprise, 
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where, naturally, any further asymptotic expansion on the adjacent saddle, i.e., for Z m (n), is to 
be evaluated precisely along the Stokes line 9 = —a nm . This discontinuity is exponentially small 
as, on the Stokes line, nWnm is real and positive. 

The Stokes discontinuity is particularly relevant to us as we later wish to identify it with 
instanton effects and, as such, we shall dwell upon it later in this section. An important thing to 
notice is that this is actually not a discontinuity of the function Z(k) but rather a discontinuity of 
the asymptotic approximation to Z(k). A trivial example is the function sinh ~. In the right half- 
plane Me z > its asymptotic behavior as z — > + is well described by | exp (-), with the error 
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being exponentially suppressed. However, should we try to rotate this asymptotic approximation 
to the left half-plane Me z < 0, it would no longer be valid as the sub dominant error — ^ exp (— -) 
will no longer be small! This is clearly a discontinuity of the asymptotic approximation we chose, 
and not of the function itself. 

We shall next see how to apply this formalism within the case of the Gamma function, which 
will later lead to a hyperasymptotic understanding of the free energies in the Gaussian and 
Penner matrix models. In particular, we shall build up our analysis in order to see how to obtain 
the multiple instanton sectors of these models straight out of the above resurgence formulae. 

5.1 Stokes Phenomena in the Gamma Function 

Consider the "Gamma partition-function" 

Z t (k)=T{k). (5.17) 
One usually defines the Gamma function via Euler's integral [80] 



f+OO 

T(k)= dww"- 1 e- w , Me(K)>0, (5.18) 
J o 

where the contour of integration is the positive real axis. The logarithm of the Gamma function, 



the "Gamma free-energy" , has a well-known representation ]80 



F r («) = logr(As) = U - ij logK-K + ilog2vr + fi(K), (5.19) 

where Q(k) is meromorphic with simple poles at K = —n, n 6 No- One then obtains asymptotic 
expansions for the Gamma function by first obtaining asymptotic expansions for the function 
f2(/c). One such familiar case is the Stirling series, which is the Poincare asymptotic expansion 

n( K ) ~ V 29 r V-r, (5.20) 

valid as \k\ — > +oo, in the sector | arg(«)| < ir. Our goal in the following is to obtain an expo- 
nentially improved version of this asymptotic expansion, in the spirit of our previous discussion 
on hyperasymptotics and Stokes phenomena, along the guidelines in pOj . 

Let us begin with the "Gamma partition-function" , applying the hyperasymptotic analysis 
in the preceding section as in pjj ]. One first changes variables as w = ne z and then re-writes 
Euler's integral fl5.18|) as 

/+oo 
-oo 

where W(z) = e z — z. The saddle points of this function are = 2irik, for k G Z, with 
W{zk) = 1 — 2-7riA;. We may now apply the machinery we previously described, with the nuance 
that one now has an infinite number of saddles. We first select the reference saddle zq = and 
define 



T q (k) = V2ttk k -2 e- K g (K), (5.22) 



Go(k) = x /—I dze-« w W- L >, Me( K )>0, (5.23) 
C (9) 
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with logt/o( K ) = ^(/c). It is not too hard [31 to identify aZZ saddles {z m } m /o as adjacent saddles 
to zo : the singulants are now 

Wo m = -2mm = 2tt \m\ exp ^(-l^+^i^Vj (5.24) 

and the steepest-descent contour from zq to each z m is that for which 6 = —aom, i-e., 

arg(e 2 - z - 1) = (-l)^^™))^ (5.25) 

As such, for m > the Stokes line is at # = £ and for m < the Stokes line is at 9 - 



2 > 



so that the imaginary axis is a Stokes line for the Gamma function. Parametrizing z = x + iy, 
this contour may also be written as cosy = (1 + x)e~ x with e~ x < if m < and greater 



than if m > (see Figure [L2|). As such, the remainder associated to the finite truncation of the 
asymptotic series for Qq{k) follows from ( |5.11 ) as 

,. n 1 +2? f f+ioo _JV-1 rO N-l } 

(5.26) 

where one should recall that the Q m (rj) are to be evaluated over the adjacent steepest-descent 
contours C m (—ao m ), e.g., for m > 0, 

g m{7]) = f± f dze -vm*)-w{z m )). (5 .27) 
V 2vr 7 Cm(f ) 

In this integral, consider the shift w = z — 2mm, where we move the contour downwards in the 
complex z-plane by 27rim. It is simple to see that the shifted contour C m (|) will now go through 
zq rather than z m and the integral becomes 



G m (v) 



ff [ dToe-^H = G (v)- (5-28) 



The exact same reasoning applies if m < 0, in which case one just shifts the adjacent contours 
upwards. Thus 

.... 1 +°° ( f+ioo N-l r0 N-l } 

4 } (*) = o— - E / dr ? e2 "" m - / . dr ? e ~ 27ri,?m ^ • ( 5 - 29 ) 

m= l wO ft J— ioo K J 

One may now address Stokes phenomena for the "Gamma partition-function" and simply 
confirm that indeed the imaginary axis in the complex K-plane is a Stokes line for Qq{k): as 
k becomes purely imaginary, one of the two integrals above will have a pole. As we have seen 
before, this leads to the discontinuities 

A ^(% =± , =±E e ±2 ^ m So(K) = T^p^So(K). (5.30) 

2 m=l 

Notice that the discontinuities are evaluated on the Stokes lines, k = ±i|/c|, and are thus al- 
ways exponentially suppressed. We will interpret these terms as nonperturbative "instanton" 
contributions to the "Gamma partition-function" . 
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Figure 12: Paths of steepest descent for the Gamma function, in the (x,y) plane, starting at zq = 0. 



As we turn to the "Gamma free-energy", the same nonperturbative corrections can be 
obtained very easily from the reflection formula 



r(«)r(-«) 



7T 



K Sin TTK 



(5.31) 



Indeed, one may now obtain the final result from a two line calculation [82]. The idea is to 
directly find the (exponentially) improved version of the asymptotic expansion ( 5.20| ), in the 
sector ^ < 9 < ir, and which we name $7 + (k). Upon the analytic continuation —k = e~ m n to 
the sector | < 9 < ir, where — § < arg (— k) < 0, one may use the expansion ( |5.19 ) for T( 
(5.31) to obtain, upon taking the logarithm of ( |5.31| ), 



-K) m 



- log (1 - e 2mK ) , 



(5.32) 



where the plus subscript refers to the sector past +~. In the sector — tt < 9 < 



-f one still 



obtains the above expression ( |5.32 ), but with e 2mK replaced by e~ 2mK instead (and the plus 
subscript naturally gets replaced by a minus subscript). Thus, in both sectors the corrections 
to the expansion ( 5.1 9| ) are always exponentially suppressed, with the discontinuity across the 
Stokes lines 9 = ±| given by 



log 1 



) = E 



+°° e ±27riKm 



m=l 



(5.33) 



which is, as expected, essentially related to the logarithm of the Stokes discontinuity for the 



'Gamma partition-function" , (5.30). It is rather tempting to understand these terms as instanton 



contributions, with "instanton action" S- 



(m) 



inst 



w( Zm ) - w(zo) = m, 



-2ix\m. 



identification to be valid, one expects that the instanton(s) with least action, ^ 



For this 
27ri and 
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= — 2m, will yield the leading contributions controlling the large-order behavior of the 
perturbative expansion ( |5.20| ), as we have discussed before in section 2.1. Setting as usual 



i + °° R 



and 



the standard large-order analysis yields 



-n<V(K) -i^e-^J^ 9 ^, (5.35) 



9=0 



o(o) „ r (gg + 6) / {1) O* 1 ^ 



(5.36) 



and, consequently, 



o (0) i /r 

9+1 - 1 +0(-\. (5.37) 



From the exact expression for the Op coefficients it immediately follows 

«L = ^l(,J + fI)) = ' +0 (i), (5.38) 

so that A = ±2-7ri just as expected from the instanton analysis. Thus, we see that indeed the one- 
instanton contributions control the large-order behavior of the perturbation theory, supporting 
the important identification of Stokes discontinuities with instanton contributions. 

An interesting property of our instanton actions is that S^" 1 ) = mS^ so that the action of 
the m-th one-instanton equals the m— th multi-instanton action. In this case, one could interpret 
the result (5.33) as being exact, in the sense of including all multi-instanton corrections, and to all 
loop orders. For this identification to be valid, one expects to fully reconstruct the perturbative 
coefficients of the "Gamma free-energy" out of its complete multi-instanton series, as in (see, 

e.g., 0) 

«+ioo 



+°° -i / p+oo rO \ 

Y^^yJ dxx 2 °- 2 e- 2 ™ x + J dxx 2 °- 2 e 2 ™ x ) . (5.39) 



-l)g +1 
2vr 



This is actually just one of the few examples on resurgent relations we have obtained in our 
earlier hyper asymptotic discussion. It is simple to see that the precise result is obtained, in full 
accordance with our multi-instanton expectations: 



(0) _ 2(-l)^(2 g -2)! ^ J_ 

m=l 



(27r)" y m 



where one just needs to use the representation of even Bernoulli numbers in terms of the zeta 
function, that we have used several times before, in order to check the result. 
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We have thus seen very clearly, at the familiar free-energy level, that the contributions 
arising from the Stokes line discontinuities are precisely the usual nonperturbative instanton 
contributions. The exact matching within the matrix model examples we consider in this work 
will be made complete in the next section. 

5.2 Instantons as Stokes Phenomena in Matrix Models 

In the previous section we have identified the Stokes discontinuities of the Gamma function 
with instanton effects of the corresponding one-dimensional integral, defining either a partition 
function or a free energy. We shall now see how this analysis carries through to the matrix models 
we are interested in. Let us start by considering the Gaussian and Penner matrix models, and 
how their instantonic sectors follow from a Stokes analysis of the corresponding free energies. 

As computed at an earlier stage, the exact free energy for the Gaussian matrix model is 
given by 

F G = ^N 2 log g s - iiV log 2tt + log G 2 (N + l), (5.41) 

and the exact free energy for the Penner matrix model (normalized against the Gaussian measure) 
is 

J=P = — log g s + ]-N 2 log g s - l-N log 2ir + log G 2 ( N + — + 1 ) - log G 2 (— + 1) . (5.42) 
9s 2 2 V g s ) \g s ) 

It is rather evident from these expressions that both cases have their asymptotic expansions 
associated to Poincare asymptotic expansions of the Barnes function. As such, the instanton 
sectors of these two matrix models will be dictated by the Stokes structure of the logarithm of 
the Barnes function. But due to the integral representation [jftj 

1 i r N 

log G 2 (N + 1) = -N log2TT- -N(N -l)+NlogT(N)- dnlogT(n) (5.43) 

2 2 J 

the Stokes structure of the Barnes free-energy is given in terms of the Stokes structure of the 
Gamma free-energy, which we have previously analyzed in greater detail. In particular 



i-N 

AlogG 2 (A + l) = JVAO(Af) -/ dnAn(n 



=±- 

^ 2 



m=l 



At the Stokes lines = ±i|A| and the discontinuities are exponentially suppressed as expected. 
It immediately follows 18 

. „ i / 2lTt S \ 2ntm 



18 Notice that the factor ^ is a genus one artifact and we drop it in the following. Indeed, when computing 
the Borel transform of the Gaussian free energy (equivalently, of the logarithm of the Barnes function), we start 
the sum at genus g = 2 in order to avoid the problematic logarithmic terms at genus zero and one, which do not 
contribute to the large-order behavior in any case . If o ne were to — incorrectly — start the sum at genus g = 1, 



while still making use of the Bernoulli expression (3.19) in order to compute this genus g — 1 term in the sum, 
one would precisely find this spurious -p? contribution. 
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and 



A/ P = -— V — ^ + ^ e § s -7— V — + ^ e 5.46 

m=l m=l x 7 

These results, where we have used t = g s N and restricted to the Stokes line at 9 = +§, and 
upon the identification of the A discontinuity with —Disc, precisely match our results for the 
full instanton sector of both Gaussian and Penner matrix models, obtained earlier in section || 
via Borel summation (or to be later analyzed via the use of trans-series methods) . We have not 
addressed the case of the Chern-Simons matrix model, as its free energy is given by the logarithm 
of the quantum Barnes function, for which we do not know of any appropriate hyperasymptotic 
framework which would allow for a derivation of its Stokes discontinuities. However, we believe 
it should be possible to study the hyperasymptotics of the quantum Barnes function in a similar 
fashion to the one above (see also ]83[|). 

Some interesting lessons may be drawn from our Stokes analysis. Because the appearance of 
the nonperturbative ambiguity of the matrix models' free energies is related to Stokes phenom- 
ena, intimately associated to discontinuities of the asymptotic approximation, we see that this 
nonperturbative ambiguity is in fact an artifact of the semiclassical, large N analysis. Clearly, the 
exact free energies, ( |5.41 ) and (5.42), are given by the logarithm of entire functions in the com- 



plex plane, with no discontinuities. The Stokes lines, and thus the instanton corrections, appear 
only at the very moment we select a particular saddle and semiclassically evaluate the partition 
functions or free energies. On the other hand, it is also clear that it is just in this semiclassical 
limit that the notion of target space in the holographically dual theory emerges. For instance, 
we have seen in an earlier section that, in the c = 1 case, the geometry of the target space arises 
from the matrix model spectral curve, which gets identified with the derivative of the FZZT 
disk partition function. Even more manifestly, in the case of the Chern-Simons matrix model 
the spectral curve of the matrix model coincides with the mirror curve of the mirror CY to the 



resolved conifold. One is led to conclude that if one is to consider the exact free energies, ( 5.41 ) 
and ( |5.42| ), as the nonperturbative definitions of the holographically dual models, then it appears 
the "exact quantum" target spaces are very different from the semiclassical ones; basically at 
the nonperturbative level the notion of target space as a smooth geometry is lost. Interestingly 
enough, this discrepancy between "semiclassical" and "exact quantum" target spaces has also 
been advocated in [j§4j[ , focusing on the example of non-critical strings with c < 1. In particular, 
Stokes phenomena was also identified therein as a source of instanton corrections. 



As we shall review in the upcoming section 6.1, in the context of matrix models the nonper 



turbative partition function is obtained by summing over all saddles of the matrix integral. In 
particular, it is the averaging over all possible semiclassical geometries that leads to the back- 
ground independence of the nonperturbative partition function, as described in [35, pLSf] - 



5.3 Smoothing the Nonperturbative Ambiguity 

As we have explained, the Stokes discontinuities are discontinuities of the asymptotic expansions 
and not of the functions under approximation, and we have made this clear in some examples. 
It so happens that in most cases one does not know the function we wish to approximate and, 
as such, one needs to devise methods to smooth the Stokes discontinuity p6| . The need for the 
smoothing should be clear: it maintains the validity of the asymptotic expansions even as we 
cross the Stokes lines. This is what we shall describe now: the universal smoothing proposed in 
[86], given by an uniform approximation involving the error function, and universally describing 
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Stokes phenomena, which naturally makes the strength of the subdominant contribution grow 
smoothly from to 1 across the Stokes line (upon where it equals |). 

The idea goes as follows. Hoping to maintain the validity of our asymptotic expansions as 
we cross Stokes lines, let us write these as 

+00 j- / \ 

„ H y 9—a nm 
g=0 m 

where we have implicitly included the remainder associated to optimal truncation within the 
infinite sum, and where we have introduced the Stokes multiplier function S m (K) weighting the 
subdominant exponentials and which will smooth the transitions across the Stokes lines. The true 
power of the Stokes multiplier function arises from the fact that, if the appropriate variables are 
chosen to cross the Stokes lines, then this function is universal, within a wide class of problems 19 . 
This appropriate "universal" variable involves the singulant W nm , specifying contours through 
adjacent saddles, i.e., specifying the location of the Stokes lines. The Stokes multiplier function 
may thus be written as |86|| 

S m (n) = i(l + erf (s nm ( K ))), (5.48) 
with erf(x) the error function and where we defined the Stokes variable 

ImfcVlWL (549) 

Observe that this is not a solution to the nonperturbative ambiguity problem as there is a choice 
of a real constant implicit in this result: the choice that before the Stokes line the exponentially 
suppressed contributions actually vanish. This is of course related to a choice of integration 
contour in the inverse Borel transform used in the calculation of the Stokes multiplier 8C ] . What 



the Stokes multiplier function does is to reduce the nonperturbative ambiguity to the choice of 
the real constant describing the intensity of the exponentially suppressed terms before the Stokes 
line, describing the crossing in an universal fashion. For a recent discussion of the nonperturbative 
ambiguity and choice of inverse Borel transform, within the matrix model context, see [14|. 

For the "Gamma partition-function" with singulant KWo m = —2ir\Km, the Stokes variable 
becomes 

, . Im (— 27riKm) /—. — - — - cos 9 n 
s m (k) = , K == = -v^MH ^=e T(1 ~ sgn(m)) 5.50 

^2 Me (-2vriKm) Vsti^d 

and the smoothing becomes implemented by the Stokes multiplier 

S m> o(K) = ^ ^1 + erf (-y/n\K\m -^= j^ , (5.51) 

where we have explicitly written the case where one crosses the upper imaginary axis. Notice 
that in the immediate vicinity of the Stokes line one has — -^== = (9 — ~ ) + • • • , simplifying 
the argument of the error function in that region. At the level of the "Gamma free-energy", 
completely analogous to the partition function analysis as the singulants are precisely the same, 



this result was interpreted in [82] as a distinct — but universal — Stokes smoothing for each ra- 
th small exponential. One may, in this light, separately understand the appearance of each 
exponential. 



3 But see |87J for a counter example. 
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For both Gaussian and Penner models, the discontinuity's singulants are the same as for the 
Gamma function and, akin to the previous discussion, one is in the presence of infinitely many 
smoothings [^|] with Stokes multipliers 

5 m >o(iV) = i(l + erf (^-|)7^V>)), (5.52) 

near the Stokes line 9 = ^. This explains the appearance of each suppressed exponential, or each 
distinct instanton contribution, in a separate but universal and smooth manner. In this way, one 
may readily obtain formal expressions for the Stokes smoothing, also yielding formal expressions 
for the free energies of these models which allow, for instance, one to cross the Stokes lines and 
reach any nonperturbative point in the complex iV-plane. 

6. Semiclassical Interpretation of Instantons 

At this stage we have a very good understanding of nonperturbative phenomena in c = 1 matrix 
models and topological strings, with all the information we have gathered both from Borel and 
Stokes analysis. However, and because at the end of the day we are analyzing large N matrix 
models in a saddle-point approximation, we would like to understand these nonperturbative 
instanton corrections directly from a semiclassical large N point of view, i.e., directly in the 
matrix model language. This is what we shall do in this section, as we provide a semiclassical 
interpretation of instantons in terms of eigenvalue tunneling, across a multi-sheeted effective 
potential, and we also suggest a spacetime interpretation for the nonperturbative effects we have 
just obtained in the preceding sections. 

6.1 Instantons as Eigenvalue Tunneling 

We shall begin by recalling how instanton effects, which control the large-order behavior of 
the 1/N expansion, are interpreted as an eigenvalue tunneling effect, within the matrix model 
context. Recalling our discussion on matrix models, in diagonal gauge the one-matrix model 
partition functions is 

Z = ^/n(^)^(A)e-i^-^). (6.1) 

As it stands this definition cannot be complete, as the integration contours for each eigenvalue Aj 
still need to be specified. It may happen that the above integral is not well-defined as a conver- 
gent, real integral, in which case the model needs to be properly defined by analytic continuation, 
i.e., by the choice of an appropriate contour in the complex plane such that the integral becomes 
convergent. Indeed, in general, the various phases of matrix models are separated by singular 
domains (in the space of complex potentials), where no large iV limit exists |2|, ||, |88|. Removal 
of the "divergent regions" where MeV(z) — > — oo as \z\ —> oo corresponds to holes in the complex 
plane, in which case one is led to decompose a generic integration path 7 on a homological basis 
of paths {71, . . . , 7 S } as 20 

s 

7 = ^Cfc7fc, (6.2) 

k=l 

20 For a polynomial potential of degree d there will be d holes in the complex plane, in which case the dimension 
of the homological basis will be s = d — 1. 
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where we shall place iVj eigenvalues on the path 7$, with arbitrary distribution {Ni} but such 
that ^2 & i= i Ni = N. As we shall make clear in the following, the coefficients Cfc may be regarded as 
theta-parameters leading to different theta-vacua ||. We may then define, with the appropriate 
symmetrizations, 



Wi „ „ N, 



€i = jj-, with i = 1, . . . , s, the free energy log Z(N\, . . . , N s ) may have a perturbative large N 
expansion 21 



For particular choices of the integration contours and particular choices of the filling fractions 

N ' 

a 21 

+00 

log Z(Ni,...,N a ) = F(N 1 ,..., N a ) = Y,N 2 - 29 F g {t), (6.4) 

where one finds the usual large-order behavior of F g ~ (2g)! rendering the topological 1/JV 
expansion asymptotic. Nonperturbative effects associated to singularities in the complex Borel 
plane are interpreted as instanton configurations as follows. 

Regarding Z(N\, . . . , N s ) as the partition function associated to a specific topological sector, 
characterized by the filling fraction e%, . . . , e s , it becomes natural to consider the general partition 
function where one sums over all possible ways of distributing the N eigenvalues 

Z(( 1 ,...Xs)= £ Ci 1 ---^'Z(N 1 ,...,N 8 ), (6.5) 

521=1 N i= N 

and where it is now clear that the Ck play the role of theta-parameters. This expression has 
been proposed by [14|, |l^] to provide a nonperturbative partition function for the matrix 
model. That ( |6.5[ ) realizes such a nonperturbative completion of the theory is made clear by 
understanding how it encodes all possible multi-instanton corrections. It was pointed out in |I5| | 
that if one is to consider the partition functions associated to two distinct topological sectors, 
with distinct fillings, {iVj} and {iV/}, one finds |T5|| 

» T »w\ dFn 

U = g s N h (6.6) 



Z(N U ...,N S ) 



implying that once one selects a reference background, {Ni} say, all other sectors are different 
instanton sectors of the matrix model. Let us consider the one-cut cubic matrix model in the 
following, in order to be a bit more concrete. The cubic potential V(z) = ^z 2 + |z 3 has two 
critical points; the maximum located at z = and the metastable minimum z = 0, as illustrated 



in Figure [L3|. There are two steepest-descent paths naturally associated to these critical points, 
70 through z = and 71 through z = z*. While the lowest energy configuration is associated 
to having all eigenvalues integrated along 70, this is an unstable configuration due to tunneling 
mediated by instanton configurations, which correspond to the integration of eigenvalues along 



71 [13]. In particular, the partition function in the one-instanton sector is given by [13] 



4 ] = ^ 41 / dx <det(*l - M'f)^ e"^. (6.7) 



21 While in general there is no topological large N expansion, there are of course some cases where this may be 
achieved: for instance in the case of degree d polynomial potential one may choose as homological basis the d — 1 
steepest-descents paths which go through each of the d — 1 critical points of the potential M, til . 
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-15 -10 -5 5 

Figure 13: The real part of the holomorphic effective potential for the one-cut cubic matrix model, in 
blue, and the simple cubic potential, in purple, for the values g = 0.1 and t = 3. 

Let us explain this expression. We removed one out of the N eigenvalues in the cut, x, and we 
are integrating it over the non-trivial saddle-point, 71. The remaining N — 1 eigenvalues are, of 
course, still integrated over the leading saddle, associated to 70, and is the zero-instanton 
partition function evaluated around this standard saddle-point. Finally, M' is an (N—l) x (N—l) 
hermitian matrix, all of its eigenvalues still integrated around the standard saddle-point in the 
zero-instanton correlation function. We refer the reader to for the details on the explicit 
computation of the quantity above. As it turns out this expression implies that, at leading order, 
the one-instanton contribution to the free energy is given by 

- fw - h w / dx (det{xl ~ M ' )2> ~-> e ~" Vi " < M > 

Z N Zj n Jxejx 

where the instanton action is fl3[| 

dzy(z), (6.9) 

with b the endpoint of the single cut C = [a,b]. This formula has an obvious semiclassical 
interpretation, as the instanton action ( |6.9[ ) is nothing but the height of the potential barrier 
under which instantons are tunneling. Furthermore, a configuration where N\ eigenvalues are 
integrated along the contour 71, and iVo = N — N\ along 70, can be naturally regarded either 
as a two-cut solution with filling fractions Nq and N±, or as a N\ -instanton excitation above 
the reference one-cut solution fll5| . It has thus become clear that the general partition function 
(|6.5| ) provides the nonperturbative completion of the matrix model since, by summing over all 
the filling fractions, it naturally encompasses all the multi-instanton configurations of the theory. 

What we shall see in the following is that, for our class of c = 1 matrix models and topological 
strings, one also has to allow for generalized integration contours which contain copies of the 
leading saddle-point and which find themselves going through different sheets of the multi-valued 
holomorphic effective potential. 
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Figure 14: The real part of holomorphic effective potentials for Gaussian, Penner and Chern-Simons 
matrix models, from left to right, respectively, when t = 0.2. 



6.2 Instantons and the Mult i Sheeted Effective Potential 

Let us now turn to the c = 1 matrix models we are considering in this work. A quick glance 
at the real part of the holomorphic effective potentials for the three models, plotted in Figure 



14, is enough to realize that there are no critical points of these potentials outside their single 
cuts. While this may not seem surprising for real and positive g s , as the asymptotic expansions 
are Borel summable and there are no nonperturbative ambiguities in the reconstruction of the 
partition functions, we also know that when allowing for imaginary values of the string coupling 
our topological expansions become non-Borel summable, with instantons controlling large order. 
As such, and given that the discussion in the previous section cannot be applied, at least not in 
a straightforward fashion, in what follows we shall have to define the instanton actions for our 
models by exploiting the structure of the holomorphic effective potentials in the complex plane. 
Let us have a closer look at the holomorphic effective potentials for Gaussian, ( 3.23| ), Penner, 



(3.38), and Chern-Simons, ( |3.52 ), matrix models. Due to the presence of either logarithmic 



or dilogarithmic functions these holomorphic effective potentials have a multi-sheeted branch 
structure in all examples, a clear feature from the structure of their Stokes lines (recall Figures 
g|| and 0). As one takes the derivative of the effective potential, in order to obtain the spectral 
curve, all logarithmic or dilogarithmic sheets artificially collapse on top of each other. As such, 
the spectral curve cannot explicitly see this structure. One way to lift this artificial degeneracy 
and make the multi-sheeted structure of the effective potentials manifest, is to repeat what we 



did in section 5.1, in the context of hyper asymptotic analysis. 

Let us start with the Gaussian model. Akin to what we did for the "Gamma partition- 
function" it is useful to first change variables from z to 



e u = z + Vz 2 - At, (6.10) 
in which case the Gaussian holomorphic effective potential gets written as 

V^ eS (u) = ^e 2u - 2t 2 e~ 2u -2tu + 2tlog2\/i. (6.11) 

In this new variable the critical points (or saddles) of the potential are located at 22 



e 



2u -. At u { ± ] = log(±2>/£) + mn, n 6 Z, (6.12) 



22 We use the standard definition of the logarithm, logz = log \ z\ + iargz, and make explicit its multi-sheeted 
structure by defining argz up to 2iri, i.e., argz| n = arg z + 2n'm with n£Z. 
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which is just the solution of a, b = e". The interesting property of the u variable is that it lifts 
the aforementioned artificial degeneracy where all sheets have collapsed on top of each other, 
and insures that one identifies all possible critical points of the holomorphic effective potential 
in its fully unfolded multi-sheet domain. One furthermore selects n = for the reference saddle 
and identifies all saddles m/Oas adjacent saddles to n = 0. 

In this new variable the multi-sheeted structure of the holomorphic effective potential is 
much simpler to visualize: one finds an infinite number of copies of the original single cut (with 
each distinct degenerate sheet parametrized by n), with endpoints u_ and u+ , and where 
^hfLf(ii) is constant on each replica of the cut. This, of course, has dramatic implications: as 
we proceed to evaluate the instanton action via A = Vh ;e ff(aio) — Vh;es(b) the sheets are no 
longer degenerate and we shall find non-zero values every time we place an eigenvalue in another 
sheet. The placing of the eigenvalue is simple and totally analogous to the previous discussion of 
eigenvalue tunneling: one removes the single eigenvalue from the endpoint of the cut and places 
it at the starting point of the "next" cut, in the following sheet (in such a way that the spectral 
curve cannot "see" any difference in the configuration, only the holomorphic effective potential 
can), an idea first suggested in [26|. We conclude that the multi-instanton action (the singulant, 
in the hyperasymptotic language we used earlier, e.g., in section ^J) is given by the difference 
between the holomorphic effective potential evaluated on the principal sheet (corresponding to 
the choice n = and denoted with a * in the following) and its value on a generic sheet, 

^ = Vgaivfi) - V£ eS (u%) = -2vrini = \ <f dzy(z). (6.13) 



In 

In the last equality we compute the instanton action as the integral of y{z) dz along j n , a non- 
contractible contour encircling the eigenvalue cut n times 23 . Analogously, the instanton action 
may be regarded as half the shift due to the additive monodromy of the logarithm present in 
Vfr eS , x \ (recall equation (|3.23|) ). The action fl6.13| ) we have just obtained from a semiclassical 
viewpoint coincides with the one previously obtained with either Borel or Stokes analysis. It is 
also important to notice that the definition we are now suggesting for the instanton action is in 
perfect analogy with the one that we have previously described when addressing the semiclassical 



derivation of the Schwinger pair-production effect, in section |2.2| . With an appropriate identi 



fication of the parameters, one may identify the semiclassical Schwinger effective action (2.25) 



with the matrix model effective potential ( |3.23 ). As such, in here one may still interpret the 



instantons as an eigenvalue tunneling process, with the sub-barrier motion consisting of multiple 
windings in the complex plane around the eigenvalue cut, as in the discussion of section 

As we move to the Penner model the holomorphic effective potential becomes more intricate, 



as displayed in ( 3.38 ), and we shall focus on obtaining the instanton action straight out of the 



monodromy shift of (3.38) around the cut C = [a, b]. There are now contributions from more than 
one logarithmic branch cut and one needs to properly consider them all. In this case, there are 
two possible ways to choose the contour: one may encircle the eigenvalue cut without crossing the 



logz cut in ( [3.38 ), in which case one is forced to cross the other two logarithmic branch cuts; or 



one may encircle the eigenvalue cut crossing the log z cut. In the first option, clockwise winding 
n times around the cut C = [a, b] one will cross two logarithmic branch cuts and obtain the global 



23 Here, 71 is the contour used to define the filling fraction in matrix models, i.e., the A-cycle. The factor of | 
insures that the action describes moving an eigenvalue from the endpoint of the cut to the beginning of the "next" 
cut in the multi-sheeted structure, i.e., that we keep the standard eigenvalue tunneling picture. 
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shift — (2t + l)27rin — 27rin = —Anin (t + l) 24 . This immediately yield the multi-instanton action 

A v n = -2vri (t + l)n=-<f dz y{z), (6.14) 

a result which agrees with the one we have previously obtained via Borel and Stokes analysis. If, 
instead, one chooses the second option, which involves crossing the logz cut, this will naturally 
yield an additional contribution of 27rim (a result we expected to find from Borel or Stokes 
analysis in any case), but it is not entirely clear to us how this relates to a cycle of the curve. 
Finally, we need to confirm the validity of our proposal within the Chern-Simons matrix 



model. In this case, the holomorphic effective potential ( 3.52 ) is rather complex, involving both 
logarithmic and dilogarithmic branch cuts. As such, we shall restrain to computing the instanton 
action from the monodromy around the cut C = [a, b]. If the cut C = [a, b] is placed on the positive 
real axis, we may choose our contour j n such that Re (z) > when contouring the cut, in which 
case one also finds IRe (£) > and both log z and log ^ will have no monodromy. We will then 
be left with the contributions of log (1 — £) and log (l — e~*£). Since both e~*£ and £ will be 
bigger than one on the right hand side of the cut, on the real axis, then both these logarithms 
will produce the same shift of 2-7ri, at each winding. Turning to the dilogarithmic dependence, 
we find a more complicated monodromy structure and we refer the reader to, e.g., ]B9[ |90[ | for a 
more in depth analysis, or to appendix |A| for a brief review of this topic. The dilogarithm Li 2 (z) 
has a branch cut starting at z = 1, and we choose it to be on the positive real axis. However, 
as one moves from the principal sheet to a generic sheet, by crossing the principal branch cut, a 
new branch cut will emerge from z = to the right. As such, we shall need two distinct integers 
in order to specify the value of Li 2 (z) on a generic sheet, in terms of its value on the principal 
sheet, Li^z). This turns out to be [^9[ |90fl 

Li 2 (z) = U%(z) + 2vrmlog(z) - Air 2 kn, n G N+, k G Z. (6.15) 

In here, the integer n counts how many times we wound clockwise in the complex plane, by 
crossing the principal branch cut, while the integer k counts how many times we crossed the 
"hidden" branch cut. Putting it all together, it follows 

Hjefrlcs = ^efflcs + 4 7rin log (e~*£) - 8ir 2 nki - 47rinlog£ + 8n 2 nk 2 = 

= v h;es\ cs ~ ^itn - 8n 2 n (h - k 2 ) . (6.16) 

Setting k = k\ — k 2 , with k\ and k 2 being the two, a priori different, windings around the 
"hidden" branch cut of the two dilogarithms in (p.52|), one may write 



dzy(z) = -4vritn - 8vr 2 nA; (6.17) 

where is a contour winding around the cut C until it reaches the (n, k) sheet. Again, on each 
(n, k) sheet we find a copy of C, with the effective potential being constant on each replica of the 
cut. The instanton action is thus given by 

A% S k = -27ritn - 47r 2 nk = - I dzy(z). (6.18) 

2 Jjn.h 



24 The reader may be puzzled by the fact the the integral of the curve around the cut, i.e., the A-cycle § A dzy(z), 
yields 2(t + 1) instead of 2t (as one would have expected from the normalization of the genus-zero resolvent at 
infinity, wo(z) ~ \ + ■ ■ ■ as z — > oo). The reason for this is that, when deforming the contour back from infinity 
to the cut C, we still have to pick the residues at z — 0. This unusual feature is, of course, due to the fact that we 
are not dealing with a polynomial potential. 
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Once again the instanton action we have obtained with our proposed semiclassical reasoning 
precisely matches the one we previously obtained from Borel or Stokes analysis. 

It is interesting to note that in all cases above the instanton action is essentially given by the 
integral of the spectral curve along the j4-cycle of the single cut, while in the cases considered in 



[13, [jjl the instanton action was given by the integral of the spectral curve along the S-cycle of 
the cut(s). This complete our matrix model derivation of the instanton action. In the following 
we focus on the spacetime interpretation of the nonperturbative effects we have discussed do far. 

6.3 Spacetime D Instanton Interpretation 

Nonperturbative effects in hermitian matrix models have been realized in terms of D-brane 
instanton effects in the holographically dual minimal models with c < 1. In particular, in || 
matrix model instantons, as described by eigenvalue tunneling, have been shown to match the 
disk contribution of ZZ branes. An analogous situation takes place in the case of c = 1 where 
the Zfk'x) disk contribution ( |3.75 ) coincides with the instanton action obtained from the MQM 



nonperturbative integral formula |43| , |25| , p9| . In this case the singularities of the curve leading 
to the ZZ brane disk amplitude are not pinched cycles, but instead are points where the curve 
has self-intersections, as described in []69| . 

As usual, let us first analyze the Gaussian matrix model, or, equivalently, the c = 1 string 
at self-dual radius. We have already observed that the MQM spectral curve is just an infinite 
covering of the hyperboloid defining the Gaussian spectral curve and identified correlators in the 
Gaussian matrix model with open amplitudes with FZZT boundary conditions. The natural next 



step is to identify the instanton action ( 6.13 ) with the ZZ brane partition function in the dual 



c = 1 theory. Moreover, considering the relation ( |3.65 ) between FZZT and ZZ boundary states 



one may write the instanton action as the difference between two FZZT branes located at the 

(n) 

branch point u+, on the principal sheet, and its replica u\ , on the n sheet, 

^FZZT/„,(n)\ 7 FZZT/„ ,* \ _ 7 ZZ 1 



A ^( = ^L\= Z™ (n^ ; ) - = Z^ = -^f dzy(z) = 27m M . (6.19) 

V l 9s J 2ig s J 7n 

A very similar story goes through for the Penner model, thus completing the identification of 
the (double-scaled) Gaussian and Penner matrix model nonperturbative effects as ZZ brane 
D-instantons, in the dual c = 1 string theory. 

Let us now turn to the case of the Chern-Simons matrix model where, much as in the local 
curve backgrounds studied in [1C], the instanton action may be interpreted, in the dual large N 



description on the resolved conifold, in terms of toric branes. In this case, the disk amplitude 
with (mirror of) toric-brane boundary conditions may be written, in terms of the mirror curve 25 , 
as BP 

Af\x) = j dsy(s). (6.20) 

It thus follows that one may write the instanton action as the difference of two disk amplitudes; 
simply set 

1 iJs / ys ys ^ys J~f nik Us 

(6.21) 



3 For backgrounds with a matrix model dual, this mirror curve precisely coincides with the spectral curve. 
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where the toric branes have been place at the endpoint u+ and at its copy, u+ ,k \ on the sheet 
labelled by n and k (recall our previous discussion). The instanton action is in this way given 
by the tension of the domain wall interpolating in between the two branes (see as well [ Q | ) . 

At this point, two comments are in order. First, when we discussed the c = 1 double-scaling 
limit we also checked that the spectral curve — alongside with all the open correlators — of the 
Chern-Simons model collapse, at the conifold point, to their corresponding c = 1 values. In 
particular, this implies that in this limit also the instanton action reduces to the c = 1 one 



which indicates, as also observed in 1C], that in the double-scaling limit nonperturbative effects 
due to toric branes reduce to the ones due to Liouville branes in minimal models. The second 
observation we want to make is that the instanton effects leading to the above action (6.21) 



involve vector multiplet moduli. This can be realized by simply noting that the instanton action 
(6.21) is given by the integral of the one-form ydz on the contour 7, which are, respectively, 



reductions of the complex structure three-form Q and of the three-cycle F for a local CY, 

A %Sk) = i J r n > r = c;-ct' fc) , (6.22) 

where the two-cycles C* and are given by line-bundles, with base points and on 

the spectral curve. The instanton action may in this way be identified with the tension of the 
domain wall interpolating in-between the two branes (see as well ]l|). 

So far we provided a spacetime interpretation of nonperturbative effects in c = 1 matrix 
models and topological strings in terms of Liouville and toric branes. However, there is a further, 
obvious, spacetime interpretation of the instanton expansion as due to BPS particle production, 
via the Schwinger effect that we have discussed earlier in this paper. Furthermore, it is possible 
to relate these particle production effects to spacetime D-brane instantons via a compactification 



to three dimensions, followed by a T-duality [91]. Recently, in |92], this connection has been 



exploited in order to study the continuity, across walls of marginal stability, of nonperturbative 
effects in type II compactifications. 



7. Trans-Series and the Toda Equation 

Another approach to the calculation of instanton corrections to the free energy of matrix models 
and topological strings, and first developed away from criticality in [jEj], deals with trans-series 
solutions in the orthogonal polynomial framework, rather than the spectral geometry as in the 
preceding section. In this section we wish to learn what this approach has to say on what 
considers Gaussian, Penner and Chern-Simons matrix models; in particular, we wish to confirm 
our instanton results in yet a novel setting. Let us begin by briefly reviewing the ideas behind 
the trans-series approach, and then apply it within our interests. 

Let us first recall that in the orthogonal polynomial formalism one may compute the partition 
function via ( 3.11| ) once one knows the recursion coefficients {r n }. Making use of the definition 



( 3.11 ) it is not too hard to obtain 

^2 = r N- (7-1) 

N 

In the continuum N — > 00 limit the coefficients r n become a function r n — > R(x, g s ) of the variable 
x = tjr G [0, t], where the function R(x, g s ) can be determined by solving the so-called pre-string 
equation, a finite difference equation obtained from the continuum limit of the recursion for the 
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coefficients {r n } Analogously, the continuum limit of equation ( |7,1| ) produces the following 
Toda-like equation for the free energy, Fn = log Zn — > F(t,g s ), 

exp (F(t + 5s , g a ) - 2F(t, g s ) + F(t - g s ,g s )) = R(t,g s ). (7.2) 

Given a solution to the recursion coefficients R(x,g s ), this equation then determines the free 
energy of our model as a solution to the Toda hierarchy. 

Let us now consider trans-series solutions (in the sense of exponential asymptotics) to the 
equations above 26 . For that, consider the trans-series ansdtz for the recursion coefficients 

+00 

R(x,g s ) = Y,C e R {e) (x,g s ), (7.3) 
e=o 

with the zero-instanton contribution given by 



hoo 

R^{x,g s ) = Y J 9l n K%l{x) (7-4) 



n=0 

and the ^-instanton contributions given by 

/ +00 



i?W(x, 5s ) = < ) (x)e-— l + ^i^x) , />!. (7.5) 



n=l 



Plugging this ansdtz into the pre-string equations one may, in principle, determine recursively 
both the instanton action A and all the loop terms Rn (x). Once this is done, one finally plugs 



the trans-series solution for R(t,g s ) on the right-hand side of equation ( |7.2| ) and solves it with 
a trans-series ansdtz for the free energy 



F{t,g s )=Y J C t F^[t,g s ), (7.6) 

where 



1=0 



+00 



F^(t,g s ) = Y J 9 2 s 9 - 2 F^ (t), (7.7) 

g=0 



and 

/ +00 



F (£) (t,<7,) = *f (t)e-i^ ^1 + £«>)J , 1. (7-8) 

Again, at least in principle, all the coefficients Fn\t) may be determined recursively. 

In [14[ this formalism has been applied and checked against a large-order analysis in several 
matrix models. Remarkably, this method appears to work beyond the context of matrix models: 
for instance, in [14], the full instanton series has been obtained for the case of Hurwitz theory, 
which is also controlled by a Toda-like equation. 



26 Notice that these trans-series solutions will only be valid in a specific region of the complex plane; as we 
change sectors and cross a Stokes line the asymptotics will change. This change will be given by a shift in the 
nonperturbative ambiguity parameter, C, as C — > C + S with S the Stokes multiplier. 
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7.1 The Trans Series Approach for c = 1 Matrix Models 

Let us now try to apply the trans-series method to solve our c = 1 models. As we have just 
learned, the first thing to do is to look for a trans-series solution to the pre-string equation. But 
it so happens that it is a common feature of all our c = 1 models that the recursion relations 
for the coefficients {r n } may be solved exactly, without a genus expansion. In particular, for 
Gaussian, Penner and Chern-Simons models we find 

r% = g s n -> R G (x,g s ) = x, (7.9) 

rl = g s n (1 + g s n) -> R p (x, g s ) = x (1 + x) , (7.10) 

r cs = q 3n (g n _ ^ with ? = ^ ^ i? cs (x, 5s ) = e 3x (e x - 1) . (7.11) 

Since we have exact solutions for the functions R(x,g s ) = R(x), without a genus expansion, there 
is clearly no asymptotics and thus no trans-series expansion for R(x). Thus, for all three cases 
that we are considering, we only have to worry about the trans-series ansdtz (78) and plug into 
the homogeneous Toda equation 

F(t + g s )-2F(t) + F(t-g s ) = 0. (7.12) 

This is a rather interesting point, also implying that all Borel poles are controlled by the Toda 
equation (and not by the pre-string equations). In hindsight this is not so surprising, as at 
criticality all our examples are in the universality class of the c = 1 string, and it is precisely 
the case that the partition function of c = 1 string theory is a r-function of the Toda hierarchy, 
satisfying the Toda equation |9«J. In some sense, the above Toda equation plays a role very 
analogous to the one played by the Painleve I equation in c = string theory, and which was 
also studied in connection to the large-order behavior of topological strings in fi"3| , |l5f , 

Let us start with the one-instanton sector. By plugging into the homogeneous Toda equation 
the I = 1 term in (78), at first order in g s one obtains 

4F«(t)e-^sinh 2 (^)=0, (7.13) 

and setting (t) ^ it follows 

sinh ( ^1) = ^ = ivrfc, keZ. (7.14) 



It is quite remarkable that the solution to the above equation already includes all fe-instanton 
actions 

A k (t) = 2iritk + a(k), (7.15) 

where we have also allowed the integration constant to depend on k. Making further use of the 
trans-series ansdtz we may compute 



+oo 



F«(t ± g s ,g s ) = Fi k \t) e"^ ( 1 + J>? ) , (7.16) 



where, for example, 



J$=lf>(t)±^(t), (7.17) 
4] = F?\t) ± £k^L( t )±^)( t ) + \^{t). (7.18) 
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We now insert these expressions into the Toda equation, and solve it perturbatively in g s . At 
second order in the string coupling it follows, 



d 2 t F[ k \t) 







F[ k \t) = $i(fc) + $ 2 (*0*- 



(7.19) 



In fact, solving the infinite chain of differential equations one perturbatively obtains, in a recursive 
fashion, one is always led to second order differential equations. All the higher order terms are 
then fixed to be 



(n+l) 



$ 2 (k)F( k \t) 



The fc-instanton contribution in Q7.8| ) is thus given by 
F^{t,g s ) = ($!(*;) + $ 2 (&)t)e * x 



n > 1. 



(7.20) 



<^ rt+1) (A;) 



n=l 



* 2 (fc) ($i(fc) + * 2 (A;)i) 



(7.21) 



It is now simple to see that the instanton series of our three models, that we have computed 
either via Borel or Stokes analysis, (|4.6|) , ( 4.16 ) and fl4.23|) , are all solutions to the homogeneous 
Toda equation 27 . This confirms, within the trans-series setting, the validity of our results. 
However, it is not equally clear how to flow in the other direction, i.e., how to obtain our results, 
(|4.6| ), (|4.16| ) and ( 1.23 ), starting from the trans-series formalism. In particular, it is not obvious 
to us how to provide enough boundary conditions in order to fix, for each distinct case, the 



integration constants of the general solution ( 7.21 ). A natural boundary condition is to impose 
matching to the c = 1 solution, in the double-scaling limit p4| . But this cannot be quite enough 
off-criticality: in fact all our distinct examples are satisfying the same difference equation with 
the same aforementioned boundary condition. Thus, one is forced to demand, as a final boundary 
condition and in order to obtain a unique solution in each case, a comparison with the large-order 
behavior of perturbation theory, in each distinct model we wish to study. So, at least for these 
models, Borel or saddle-point methods seem more powerful and efficient ways to proceed. 



7.2 A Comment on Parametric Resurgence 

In many cases it is not possible to resum the Borel transform of an asymptotic series, or even to 
explicitly compute this Borel transform, and thus one cannot locate the singularities in the Borel 
complex plane (some of which, in particular, control the large-order behavior of the perturbation 
theory). However, it may be the case that, even if we are not able to resum the Borel transform, 
we may know that the asymptotic series arises as a solution to a finite difference equation (where 
we are, of course, interested in the example of the Toda equation). In this case we may still 



obtain some interesting information, as shown in [95]. Let us quickly apply [95] to our problem. 
Consider the asymptotic series 



+oo 



F(t,g s )^J29s n Fn(t), 



(7.22) 



n=0 



27 Because the Toda equation is linear, linear combinations of (7.21) are also a solution to the problem, in 
particular linear combinations where a also depends on an integer m and we sum over m € Z as in the conifold. 
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which we take as a perturbative solution to the finite difference equation 

F(t + g s ) - 2F(t) + F(t - g s ) = g 2 s G(t) O 4sinh 2 F(t) = g 2 s G(t). (7.23) 



Then, it is possible to formally solve this differential equation [35], with the formal solution being 



expressed in terms of the power-series coefficients of the function 

2 +°° 

ffW = - 1 + I^i)=2>* 2 " +2 - (7 - 24) 

\ 2 ) n=0 

Furthermore, the poles of the Borel transform i3[.F](£) may be related to the poles of the Borel 
transform of H(x). In particular [^] this implies that in this case the Borel poles are located at 
2iTmt, with n G Z. It is indeed the case that the instanton actions for the three matrix models 
we have studied in this paper are of this kind. 

8. Conclusions and Outlook 

In this paper we have addressed the nonperturbative structure of topological strings and c = 1 
matrix models, focusing on Gaussian, Penner and Chern-Simons matrix models together with 
their holographic duals, c = 1 minimal strings and topological strings on the resolved conifold. 
Making use of either Borel or Stokes analysis, we have uncovered the nature of instanton effects 
in these models, and have further explored the relation of these nonperturbative phenomena to 
the large-order behavior of the 1/N expansion. While this builds up on previous work along the 
same direction ||i"3| , 14, l||, IS], clarifying the discussion for a big class of models, we believe there 



is still much work to be done. In particular, let us end by listing several issues raised in this 
paper which we believe deserve further and immediate investigation (in no particular order): 

• While our results were checked in the trans-series formalism, we have seen that it is not 
obvious how to obtain these results for the multi-instanton series starting straight from 
within the trans-series set up. In particular, the choice of boundary conditions is not 
completely clear. It would be interesting to further explore and better understand the 
trans-series ansatz in this context, possibly solving the questions we have just mentioned. 

• In the Gaussian and Penner models we have explicitly shown that instanton effects may be 
understood as Stokes phenomena for the logarithm of the G2(z) Barnes function. However, 
we could not say much along these lines for the Chern-Simons model, as we do not know 
of any appropriate hyper asymptotic framework for the quantum Barnes function, G q {z). 
It would be very interesting to study the hyperasymptotics of G q {z) and show that Stokes 
phenomena in this case is also related to the instanton effects of the Chern-Simons model. 

• In our matrix model derivation of the instanton effects, in terms of eigenvalue tunneling, we 
have obtained the instanton action expressed in terms of the spectral curve (a cycle of the 
y dz one-form) . It would be rather interesting to extend this calculation in order to contem- 



plate loop corrections. Indeed, in [13, 15], higher loop terms around the multi-instanton 
configurations were computed, in terms of matrix model open correlators. Extending that 
calculation to the present set up would be very interesting, since these correlators can be 
computed entirely in terms of the spectral curve and, as such, this formalism could be 
extended to other topological string scenarios where a dual matrix model description is 



- 60 - 



not available. In particular, this could allow for a direct understanding — from a spectral 
geometry point of view — of why the loop expansion around an one/multi-instanton config- 
uration truncates in all our examples. Furthermore, this would also provide an explanation 
of the instanton expansion in terms of open string amplitudes with either Liouville or toric 
boundary conditions, for the c = 1 string and the resolved conifold, respectively. 

Much of our Borel analysis was very much related to the existence of a GV integral repre- 
sentation for the free energies of our models. Since also on a generic CY background the 
topological string free energy admits a GV integral representation, one is led to wonder if 
our approach may be applied to other, more general cases. Recall that on a general CY 
threefold X, the topological string free energy is given in terms of GV integer invariants by 
the expansion 

When r = the zeros of the sine will be poles of the integrand and this contribution to the 
total sum will, as in the case of the resolved conifold which we address in the paper, yield 
a nonperturbative contribution to the free energy. However, for generic CY backgrounds, 
one will also have to consider the summation over the Kahler classes {di}, as weighted by 
n^'. In this way, one is led to write (we neglect the pole at zero) 

D isc as.) = - ' g ft g £ fQ*fm±M + 1\ e -^, +i ,„> 

ys di = l n=lm6Z V / 

(8.2) 

Higher terms with r > in (|8.1| ) have no poles in the complex plane and will only contribute 
to the nonperturbative corrections through the residues at infinity (recall (|2.8[) and our 
discussion in the appendix). What role they might play is beyond the scope of our analysis. 

It seems likely that nonperturbative corrections obtained from the GV representation ac- 
tually provide the full nonperturbative corrections to the topological string free energy, in 
those cases where the number of GV invariants rir(X) is finite (as in our example of the 



resolved conifold). Indeed, in these cases (8.2) may provide the complete tower of nonper- 



turbative corrections to the topological string free energy, as ( |8.1| ) is basically given by a 
finite sum, not a power series expansion. In particular, the number of njf invariants is 
finite for non-singular curves of any genus |)6| , for rational curves with nodal singularities 



97], and for the configurations studied in |9^[ (among which are the CY threefolds which 
are Ak-type ALE spaces, times C). Further notice that in fl8.2Q the loop expansion around 
multi-instanton configurations truncates. In this case, and if indeed ( |8.2|) turns out to be 
the full answer for backgrounds with a finite number of GV invariants, then it must also 
be the case that these backgrounds will display no non-trivial large-order behavior in their 
multi-instanton sectors. Finally, observe that the A^-type ALE backgrounds have also 



been studied in the context of the OSV conjecture [21|, in [99]. The OSV conjecture [21] 
relates the topological string partition function to the partition function of a configuration 
of branes in type II string theory, giving rise to a four-dimensional BPS black hole. The 
brane partition function was further suggested to provide a nonperturbative completion of 
topological string theory. For A^-type ALE spaces, times C, it would thus be interesting 
to compare OSV and Schwinger completions. 
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When the number of GV invariants rir(X) is infinite, it seems very unlikely that the GV 
integral representation can still provide the full set of nonperturbative corrections to the 
topological string free energy. Indeed, it is now the case that ~ e d , for large d, and it 



seems to be the case that (8.2) cannot provide the complete nonperturbative information. A 
particularly interesting example to further explore this issue is that of the local curve, which 
indeed has an infinite number of GV invariants. The good news is that this background 
has been extensively studied in [ 13 1 , with instanton configurations identified and checked 
against the large-order behavior. As such, it would be extremely interesting to analyze the 
relation between the nonperturbative corrections arising due to the GV representation, with 



the ones derived in [13[. In particular, the analysis of 13, 15] seems to indicate that, for the 
local curve, the loop expansion around multi-instanton configurations will not truncate, 
which is to say that (|8.2| ) cannot be the full correct answer. On the other hand, and as 
we have already remarked, the local curve is in the universality class of 2d gravity with 
c = 0, and as such it should better be understood as a first step to understand backgrounds 
with an infinite number of GV invariants. A second natural step in this direction would 
be looking at the case of local P 2 , a c = 1 toric geometry with an infinite number of GV 
invariants. Indeed, the free energy of this model can be computed very efficiently to all 
genus, by means of direct integration of the holomorphic anomaly equations as shown in 



100], and would thus provide for a natural testing ground for our aforementioned questions. 



Besides an explicit check against the large-order behavior of the theory, another way to test 
the Schwinger completion of topological string theory, for generic backgrounds, would be 
to study the modular properties of its nonperturbative free energy. In fact, it is expected 
that modular invariance may be recovered at the nonperturbative level fl9|| , an important 
issue also in the context of large N dualities. Clearly, for the case of the resolved conifold, 
which we studied in this paper, there are no constraints arising from modularity since the 
moduli space is trivial, with the mirror geometry having genus zero, but this will not be 
the case for, e.g., local P 2 , whose mirror geometry has a spectral curve of genus one. Notice 
that backgrounds with a finite number of GV invariants seem not to give rise to mirror 
geometries with spectral curves of genus one and, as such, modularity should become an 
issue precisely when flS.2| ) ceases to be the full correct answer. In particular, modularity 
may play a key role in order to understand exactly what type of information is ( |8.2[ ) missing 
in general backgrounds, and these issues should be addressed in future work. 



Acknowledgments 

We would like to thank Jacopo Belfi, Andrea Brini, Luca Griguolo, Jose Mourao, Nicolas Orantin, 
Christoffer Petersson, Domenico Seminara, Jorge Drumond Silva, Angel Uranga, Marcel Vonk 
and, specially, Marcos Marino, for useful discussions, comments and/or correspondence. RS 
would like to thank CERN TH-Division for hospitality, where a part of this work was conducted. 



- 62 - 



A. The Polylogarithm: Branch Points and Monodromy 



This appendix is devoted to the study of the polylogarithm, with emphasis towards its branch 
points and monodromy. This function is defined by 

71=1 

and on its principal sheet it has a branch point at z = 1 leading, by convention, to a branch 
cut discontinuity in the complex z plane running from 1 to infinity. As one starts "exploring" 
the multi-sheeted structure of the polylogarithm and moves off its principal sheet, one finds that 
there exists another branch point, at z = 0. In this case, the resulting monodromy group will 
be generated by two elements, acting on the covering space of the bouquet S 1 V S 1 of homotopy 
classes of loops in C\{0, 1}, passing around the branch points z = or z = 1. For further details, 
we refer the reader to the very thorough explanations that can be found in, e.g., |K| |3S| ]. 

For our purposes in this paper, a simple analysis in terms of explicit topological language 
will suffice. Let mi represent the homotopy class of all loops based at some point z in C, which 
wind once, clockwise around the branch point at z = 1. The action of mi on the polylogarithm 
has the effect of carrying this function from one sheet to the next. It was shown in that one 
may write 

m 1 -U s (z)=Li s (z)-A 1 , (A.2) 

where Ai is a function, whose specific form is not important at the moment, but which includes 
a logarithm with a branch point at z = 0. This implies that, after acting once with mi, one finds 
oneself on a sheet which has a branch cut discontinuity running from to minus infinity. If we 
now let mo represent the homotopy class of all loops based at some point z in C, which wind 
once, clockwise around this new branch point at z = 0, its action on the logarithm is the familiar 
one: 

mo • log z = log z + 2-7ri. (A.3) 

Now, because the principal sheet of the polylogarithm has no branch point at z = 0, it simply 
follows 

m • U s (z) = U s (z). (A.4) 
If one now winds with mi in the opposite direction, one is led to write instead 

m^ 1 • hi s (z) = Li s (z) - A_i, (A.5) 



where again A_i is a function we shall leave unspecified; see [90] for details. If mi is to be properly 
considered the group-theoretic inverse of mi it better be the case that mi • m^ 1 = 1 = m^ 1 ■ mi, 
when acting on hi s (z). This immediately implies, e.g., m^ 1 Ai = — A_i, where we recall that Ai 
includes the standard logarithmic branch cut discontinuity starting off at z = 0. This relation 
thus seems odd, as the logarithm has no branch point at z = 1 and there should be nothing to 



wind around. This is a subtle point, further explained in |X|, and it should be stressed that it is 
the joining of polylogarithmic and logarithmic cuts that causes this effect. One way to capture 
the idea of there being no obstruction for the logarithm at z = 1 is [ p0| ] to define group elements 
gi = mi ■ jjiq 1 and go = mo, such that gi ■ log z = log z and 

go • Li s (z) = H s (z), (A. 6) 

gi-Li s (z) = U s (z)-Ai. (A.7) 
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We shall explore this monodromy group for the dilogarithm in the following. 

The free combinations of powers of the two generators go and g± form the monodromy group 
of the polylogarithm. If s is a positive integer, this monodromy group has a finite-dimensional 
representation with dimension s + 1. A particularly well known case is the dilogarithm s = 2, 



also further discussed in [89]. In this case, the monodromy group is the discrete Heisenberg group 
1 90]. In particular, one finds 



A n = 2ni (log z + 27ri (n - 1)) 



(A. 



As such, repeated applications of go and g\ will only result in linear combinations of the diloga- 
rithm Li2(z), the logarithm logz, and the identity operator. Indeed, one could further take each 
of these three elements as a basis of a three-dimensional vector space, ei = 47r 2 , e2 = — 27ri log z 
and e3 = Li2(z), in which case the matrix representation of the monodromy group would become 
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10 
1 1 
1 



and 



9o 



1 1 
10 
1 



(A.9) 



These two matrices are in fact the generators of the discrete Heisenberg group %3(Z), see I9 1 
for full details on this discussion. 

In this paper we are interested in the following action (see [pC( | for any missing details) 



g% ■ gl ■ Li 2 (z) = 5o ■ (Li 2 (z) - n A x ) = 

= Li2(z) — n Afc+i = \j\ 2 {z) — 2irin logz + 4ir 2 kn, 

and we make use of this result in the main body of the paper. 



(A.10) 



B. Dispersion Relation for Topological Strings 

In this appendix we address the Cauchy dispersion relation ( |2.8| ) for the case of general topological 
string theories on a CY threefold X , whose free energy is given in terms of GV integer invariants 
by the expansion ( |2~37l ), and which we recall in here as 

+00 +00 +00 1 9r-_9 

W = ^^„W W ^I(2sin|) " e- M (B.l) 

r=0 di=l n=l 

In particular, we wish to evaluate 

I (B.2) 
J(oo) 2tti w- g 

and show that this vanishes in the case of the resolved conifold, a result we have used in the 
main body of the paper. Notice that one cannot compute the pole at infinity straight: infinity is 
an essential singularity of the integrand in the GV representation and, as such, residue calculus 
does not apply. Nonetheless, for r = one may decompose the contour Coo into the sum of two 
contours in upper and lower hemispheres, C+ and C_, respectively, 

F dw 1 1 [ dw e inw [ dw e-' mw 

JCoo 27ri w ~9 (2sin^) 2 Jc+ 27ri (to - g) (1 - e inw ) 2 Jc- 27ri (w - g) (1 - e ~ inw ) 2 ' 



- 64 - 



and then use Jordan's lemma — applicable as 



lim max 

i?-H- oo0e[O,7r] 



(Re ie -g) (l-e inRei9 Y 



(B.4) 



in the upper hemisphere, and analogously in the lower — in order to find that this implies that the 



integral vanishes at infinity, and it thus follows that for the resolved conifold, where only n 



(i) 
o 



1 

is non-zero, the ( |B.2| ) contribution indeed vanishes. For more complicated CY threefolds with 
r > 1 it seems rather likely that there will be a Cauchy contribution at infinity, and a complete 
analysis of this situation is beyond the scope of the present work. 
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